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Abstract 

This paper consists of three parts. In part I, we microscopically derive Ginzburg- 
Landau (GL) theory from BGS theory for translation-invariant systems in which mul¬ 
tiple types of superconductivity may coexist. Our motivation are unconventional su¬ 
perconductors. We allow the ground state of the effective gap operator Kt^ + V to be 
n-fold degenerate and the resulting GL theory then couples n order parameters. 

In part II, we study examples of multi-component GL theories which arise from an 
isotropic BGS theory. We study the cases of (a) pure d-wave order parameters and (b) 
mixed (s -|- d)-wave order parameters, in two and three dimensions. 

In part III, we present explicit choices of spherically symmetric interactions V which 
produce the examples in part 11. In fact, we find interactions V which produce ground 
state sectors of + V of arbitrary angular momentum, for open sets of of parameter 
values. This is in stark contrast with Schrodinger operators —-|- V, for which the 
ground state is always non-degenerate. Along the way, we prove the following fact 
about Bessel functions: At its first maximum, a half-integer Bessel function is strictly 
larger than all other half-integer Bessel functions. 
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1 Introduction 

Since its advent in 1950 [19], Ginzburg-Landau (GL) theory has become ubiquitous 
in the description of superconductors and superfluids near their critical temperature 
Tc- GL theory is a phenomenological theory that describes the superconductor on a 
macroscopic scale. Apart from being a very successful physical theory, it also has a 
rich mathematical structure which has been extensively studied, see e.g. [11, 14, 21, 40] 
and references therein. Microscopically, superconductivity arises due to an effective 
attraction betweeen electrons, causing them to condense into Cooper pairs. In 1957 
Bardeen, Cooper and Schrieffer [6], were the first to explain the origin of the attractive 
interaction in crystalline s-wave superconductors. By integrating out phonon modes, 
they arrived at their effective “BCS theory”, in which one restricts to a certain class of 
trial states now known as BCS states. In 1959, Gor’kov [20] argued how the microscopic 
BCS theory with a rank-one interaction gives rise to the macroscopic GL theory near 
Tc- An alternative argument is due to de Gennes [12]. 

The first mathematically rigorous proof that Ginzburg-Landau theory arises from 
BCS theory, on macroscopic length scales and for temperatures close to Tc, was given 
in [16] under the non-degeneracy assumption that there is only one type of super¬ 
conductivity present in the system. The derivation there allows for local interactions 
and external fields and hence applies to superfluid ultracold Fermi gases, a topic of 
considerable current interest. 

In the present paper, we use the same formalism as in [16] and study microscopically 
derived Ginzburg-Landau theories involving multiple types of superconductivity for 
systems without external fields. 

We first discuss the main result of part I, which forms the basis for the applications 
in parts II and III. Afterwards, we discuss the physical motivation for studying multi- 
component GL theories and the extent to which our model applies to realistic systems. 
The introduction closes with a description of the main results of parts II and III. 

1.1 Main result of part I 

As in [16], we employ a variational formulation of BGS theory [5, 33] with an isotropic 
electronic dispersion relation. We use previous rigorous results about this theory in the 
absence of external fields [15, 22, 24, 25]. Particularly important is the result of [22] 
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that the critical temperature Tc can be characterized by the following linear criterion. 
Tc is the unique value of T > 0 for which the “effective gap operator” 


Kt + F(x) 


tanh 


+ y(x) 


has zero as its lowest eigenvalue. Here V is the electron-electron interaction potential. 
Throughout the microscopic derivation of GL theory in [16], it is assumed that zero is 
a non-degenerate eigenvalue of Kt^ -t- V. For radially symmetric V, this means that 
the order parameter is an s-wave, i.e. it is spherically symmetric. 

The main result of part I, Theorem 2.10, is that for systems without external 
helds the microscopic derivation of GL theory also holds when the eigenvalue is degen¬ 
erate of arbitrary order n > 1. (A general argument shows that always n < oo.) The 
arising GL theory now features precisely n order parameters ipi, ... ,ipn- It turns out 
that one can use the same general strategy as in [16]. 

In fact, one can classify approximate minimizers of the BGS free energy via the GL 
theory. Given an orthonormal basis {oi,... ,an} of ker(Ar 7 ’^ -|- V), Theorem 2.10 (ii) 
says that, near the critical temperature, the Cooper pair wave function a of a BGS 
state of almost minimal free energy (i.e. the Cooper pair wave function realized by the 
physical system) is approximately given by a linear combination of the {ai,..., On} of 
the form 



where the “amplitudes” V’l > • • • > V’n £ C almost minimize the corresponding GL func¬ 
tion. 

The results of [16] allow for the presence of weak external fields which vary on 
the macroscopic scale. A key step is to establish semiclassical estimates under weak 
regularity assumptions. We emphasize that in our case the system has no external 
helds and is therefore translation-invariant. This simplihes several technical difficulties 
present in [16]. In particular, the semiclassical analysis of [16] reduces to an ordinary 
Taylor expansion. The result of the expansion is stated as Theorem 5.3 and we give 
the simplified proof for the translation-invariant situation. We do this (a) to obtain 
optimal error bounds and (b) to hopefully make the emergence of GL theory more 
transparent in our technically simpler situation. 


1.2 Physical motivation 

Background. The degenerate case corresponds to systems which have multiple or¬ 
der parameters, i.e. which can host multiple types of superconductivity. Physically, this 
situation occurs e.g. for uneonventional supereonductors. By definition, these are ma¬ 
terials in which an effective attractive interaction of electrons leads to the formation 
of Cooper pairs, but the effective attraction is not produced by the usual electron- 
phonon interactions. (Identifying the underlying mechanisms is a major open problem 
in condensed matter physics.) 
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Two important classes of unconventional superconductors are the layered cuprates 
and iron-based compounds, typically designed to have large values of Tc (“high-temperature 
superconductors”). Many of these materials possess tetragonal lattice symmetry, though 
the prominent example of YBCO has orthorhombic symmetry. There is strong exper¬ 
imental evidence for the occurrence of d-wave order parameters in these materials, in 
contrast to the pure s-wave order parameter in conventional superconductors. More 
precisely, phase-sensitive experiments with Josephson junctions [28, 29, 46, 47, 51] have 
evidenced the presence of a da. 2 _j^ 2 -wave order parameter (for tetragonal symmetry) and 
of mixed (s -|- d^ 2 _y 2 )-wave order parameters (for orthorhombic symmetry). 

There also exist proposals of d-wave superfluidity for molecules in optical lattices 
[30]. 

Multi-component Ginzburg—Landau theories. On the theoretical side, one 
of the most important tools for studying unconventional superconductors are multi- 
component Ginzburg-Landau theories [2, 7, 26, 34, 39, 44, 48, 49, 52, 53]. Many of 
these papers study the symmetry properties near the vortex cores in two-component 
GL theories. A very common example is a GL theory with (s -|- d^ 2 _y 2 )-wa,ve order 
parameters; this case has also been studied mathematically in [13, 27]. The effect of 
an an anisotropic order parameter on the upper critical field was studied in [32]. 

Another avenue where two-component GL theories have been successful is in the 
description of type 1.5 superconductors [4, 9, 42]. These are systems in which the 
magnetic field penetration depth lies in between the coherence lengths of the different 
order parameters (of course this effect only manifests itself in an external magnetic 
field). 

Microscopically derived GL theories. In many of the papers cited above, 
the GL theories that are studied are first obtained microscopically by using Gor’kov’s 
formal expansion of Green’s functions. The advantage of having a microscopically 
derived GL theory is that it has some remaining “microscopic content”. By this we 
mean: 

1. One can directly associate each macroscopic order parameters with a certain 
symmetry type of the system’s Cooper pair wave function. Therefore, if we can 
classify the minimizers of the microscopically derived GL theory, we understand 
exactly which Cooper pair wave functions a can occur in the physical system in 
configurations of almost minimal free energy. 

2. One has explicit formulae for computing the GL coefficients as integrals over 
microscopic quantities. 

The first point is expressed by (1.1) above and is therefore a corollary of Theorem 
2.10. The second point is represented by formulae (2.20),(2.21) in Theorem 2.10. 

While the papers cited above provide important insight about the vortex structure 
in unconventional superconductors, they are restricted in that the GL theories are ob¬ 
tained using the formal Gorkov procedure and that almost exclusively two-component 
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GL theories are studied. Our Theorem 2.10 provides a rigorous microscopic deriva¬ 
tion of n-component GL theories with n arbitrary starting from a BCS theory with an 
isotropic electronic dispersion. 

Physical assumptions of our model. We discuss the main physical assumptions 
of our model and the resulting limitations in its applicability to realistic systems. 

(a) Translation-invariance. We view the degenerate translation-invariant systems as 
toy models for multi-component superconductivity. We believe that the examples 
of multi-component GL theories studied in part II are already rich enough to show 
that the translation-invariant case can be interesting. From a technical perspec¬ 
tive, translation invariance yields major technical simplihcations. In particular, 
the semiclassical analysis of [16] reduces to a Taylor expansion. 

(b) BCS theory with a Fermi-Dirac normal state. There are two assumptions here: 
First, we start from a BGS theory (meaning a theory in which electrons can 
form Gooper pairs and which restricts to BGS-type trial states). The question 
whether such a theory can be used to describe unconventional superconductors is 
unresolved [34]. Second, we work with a BGS theory for which the normal state 
is given by the usual Fermi-Dirac distribution. Most realistic unconventional 
superconductors are strongly interacting systems with a non-Fermi liquid normal 
state [34, 41]. 

(c) Isotropy. We study a BGS theory in which the electrons live in the continuum 
and have an isotropic dispersion. Many of the known examples of unconventional 
superconductors are layered compounds which are effectively two-dimensional. 
When we say that their order parameter has dx 2 _y 2 -wave symmetry, then this only 
means that it has a four-lobed shape similar to that of k^ — ky for —tt < kx,ky < tt, 
but its precise dependence on kx, ky depends on the symmetry group of the two- 
dimensional lattice [34]. Order parameters of the form k"^ — ky have been studied 
as a hrst approximation to unconventional superconductors, see e.g. [39, 49, 52]. 
For the examples in part II, we consider a spherically symmetric interaction po¬ 
tential, resulting in a fully isotropic BGS theory. Gonsequently, the d-wave order 
parameters that we consider are the “usual” ones, known from atomic physics 
(see section 3.1). By isotropy, all the d-waves (there are two in two dimensions 
and five in three dimensions) are energetically equal. The examples in part II 
show that even this isotropic microscopic theory can lead to rather rich coupling 
phenomena of anisotropic macroscopic order parameters, as we discuss next. 

(d) Spin singlet order parameter. We restrict to order parameters which are singlets 
in spin space. This is indeed the case for unconventional superconductors [34], but 
it excludes systems with p-wave order parameters such as superfluid Helium-3. 

1.3 Main results of part II 

In part II, we compute the n-component GL theories that arises from the BGS theory 
according to Theorem 2.10 for several exemplary cases. For each situation, we make 
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some observations about the minimizers of the GL energy and their symmetries and 
give a physical interpretation. 

Throughout part II, V is assumed to he spherically symmetric, so the BCS theory 
becomes fully isotropic. The order parameters can then be described by the decom¬ 
position into angular momentum sectors (see section 3.1) and we consider the case of 
pure d-wave and mixed (s -|- d)-wave order parameter. Here and in the following, we 
write “GLn” for “n-component Ginzburg-Landau theory”. The dimension D will be 
either two or three. 

(i) Let D = 3. Assume the Cooper pair wave function is a linear combination of 
the five linearly independent d-waves with a given radial part. Theorem 3.1 
explicitly computes the microscopically derived GL5 energy and gives a full de¬ 
scription of all its minimizers. Surprisingly, the GL5 energy in three dimensions 
exhibits the emergent symmetry group 0(5), see Corollary 3.3 (i), which is con¬ 
siderably larger than the original 0(3) symmetry group coming from the spherical 
symmetry and reflection symmetry of V. 

(ii) Let 0 = 2. Assume the Cooper pair wave function is a linear combination of the 
two linearly independent d-waves with a given radial part. Theorem 3.5 explic¬ 
itly computes the microscopically derived GL2 energy and gives a full description 
of all its minimizers. We find that the {d,j. 2 _y 2 ,dxy) order parameter must be of 
the form (■0, TiV') with \^\ minimizing an appropriate GLl. In particular, the 
minimizers of this GL2 form a double cover of the minimizer of a GLl. 

(hi) Let 0 = 3. Assume the Cooper pair wave function is a linear combination of the 
five linearly independent d-waves with a given radial part and the s-wave with 
another given radial part. Theorem 3.7 explicitly computes the microscopically 
derived GL6 energy. It also gives a simple characterization of the parameter values 
for which the pure d-wave minimum is always unstable under s-wave perturbations 
and of the parameter values for which, vice-versa, the pure s-wave minimum is 
unstable under d-wave perturbations. As a consequence, we give parameter values 
for which s- and d-waves must couple non-trivially to be energy-minimizing. 

We also consider the mixed (s -|- d)-wave case in O = 2 dimensions. The result is 
presented in Remark 3.9 (v) for brevity. 

1.4 Main results of part III 

Recall from the discussion of part I above, that the candidate Cooper pair wave func¬ 
tions are the ground states of the effective gap operator iLj’c + R- A priori, it is not 
at all clear that the fully isotropic BCS theory can produce ground state sectors of 
Htc + R which are not spherically symmetric. In particular, it is not clear that the 
examples considered in part II actually exist. 

In fact, if Rtc is replaced by the Laplacian — we have a Schrodinger operator 
and under very general conditions on the potential R, the Perron-Frobenius theorem 
implies that the ground state is in fact non-degenerate, see e.g. Theorem 11.8 in [35]. 
For spherically symmetric R, this means the ground state is also spherically symmetric 
(“s-wave”). 


6 



In part III, we remedy this by exhibiting examples of spherically symmetry poten¬ 
tials V such that the ground state sector of + V can in fact have arbitrary angular 
momentum. These potentials will be of the form 

= -A(I(|x| - R) 

in three dimensions. Here A and R are positive parameters. The result holds for open 
intervals of the parameters values, so it is “not un-generic”. 


2 Part I: Microscopic derivation of GL theory 
in the degenerate case 


2.1 BCS theory 

We consider a gas of fermions in with 1 < D < 3 at temperature T > 0 and 
chemical potential /r G M, interacting via the two-body potential H(x). We assume 
that V(x) = I/(—x) is reflection symmetric. We do not consider external fields, so the 
system is translation-invariant. A BCS state T can then be conveniently represented 
as a 2 X 2 matrix-valued Fourier multiplier on of the form 


r(p) 


7(P) a(p) \ 
«(p) 1 - 7(P) ) ’ 


( 2 . 1 ) 


for all p G Here, 7 (p) denotes the Fourier transform of the one particle-density 
matrix and 3(p) the Fourier transform of the Cooper pair wave function. We require 
3(p) = 3(—p) and 0 < F(p) < 1 as a matrix, which is equivalent to 0 < 7 (p) < 1 
and |3(p)p < 7 (p)(l — 7(p)). The BCS free energy per unit volume reads, in suitable 
units 

j (p 2 -^)7(p)dp-r5[F] + / H(x)|a(x)pdx, (2.2) 

Jrd J]rd 

where the entropy per unit volume is given by 


5[F] 


Trc2 


r(p) logf(p) 


dp. 


(2.3) 


Remark 2.1 (BCS states), (i) In general [5, 16], 5?7(2)-invariant BCS states are rep¬ 
resented as 2 X 2 block operators 

7 a_\ 

\a 1-77 

where 7, a are operators on L^(]R^) with kernel functions 7(x, y) and q;(x, y) 
in L^(M^) © L^(M^). Since 0 < F < 1 is Hermitian, 7(x,y) = 7(y,x) and 
q;(x, y) = Q;(y,x). In the translation-invariant case considered here, these kernel 
functions are assumed to be of the form 7(x —y) and a(x —y). Since convolution 
by 7, a becomes multiplication in Fourier space, we can equivalently describe the 
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BCS state by its Fourier transform F defined in (2.1) above. In the translation- 
invariant case, the symmetries of 7 , a turn into the relations 7 (x) = 7 (—x) and 
q;(x) = a(—x) or equivalently 7 (p) = 7 (p) and S(p) = 3(—p). Finally, since we 
are interested in states with minimal free energy, we may also assume 

7(P)=7(-P) (2.4) 

and this was already used on the bottom right element in (2.1). To see this, let 
F be a BCS state not satisfying (2.4), set rr(p) := r(—p) and observe that 

jrBcs j'TtIr) < ^Bcs(r) 

by strict concavity of the entropy and reflection symmetry of all terms in . 

(ii) Note that a(x, y) = Q;(y, x) means that the Cooper pair wave function is symmet¬ 
ric in its arguments. To obtain a fermionic wave function, we would eventually 
tensor a with an antisymmetric spin singlet. Since ol is reflection-symmetric in 
the translation-invariant case, a. must be of even angular momentum if V is radial. 
The restriction to symmetric a is a consequence of assuming SU{2) invariance 
in the heuristic derivation of the BCS free energy [ 22 , 33]. This means the full 
Cooper pair wave function must be a spin singlet and so its spatial part a must 
be symmetric. Note that this excludes systems, e.g. superfluid Helium-3, which 
display a p-wave order parameter. 

(hi) For more background on the BCS functional, in particular a heuristic derivation 
from the many-body quantum Hamiltonian in which one restricts to quasi-free 
states, assumes SU{2) invariance and drops the direct and exchange terms, see 
[33] or the appendix in [22]. Recently, [ 8 ] justified the last step for translation- 
invariant systems by proving that dropping the direct and exchange terms only 
leads to a renormalization of the chemical potential p, for a class of short-ranged 
potentials. 

We make the following technical assumption on the interaction potential. 

Assumption 2.2. We either have V G LP'^(M^) with pv = 1 ior D = 1, 1 < pv < 00 
for D = 2 and pv = 3/2 for D = 3, or we have 

F(x) = Fa,r(|x|) := -A5(|x| - R), (2.5) 

when H = 1, 2, 3 and A, R > 0. 

We note 

Proposition 2.3. A potential V satisfying Assumption 2.2 is infinitesimally form- 
bounded with respect to — V^. 

We quote a result of [22], which provides the foundation for studying the variational 
problem associated with . Define 

V := {f as in (2.1) : 0 < F < 1, 7 G , (1 + p^) dp), a G R,V(M^)} 

with R/y^(M^) = {a G Ff^(M^) : a(x) = a(—x) a.e.}. 
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Proposition 2.4 (Prop. 2 in [22]). Under Assumption 2.2 on V, the BCS free energy 
(2.2) is bounded below on V and attains its minimum. 

The physical interpretation rests on the following 

Definition 2.5 (Superconductivity). The system described by is superconduct¬ 
ing (or superfluid, depending on the context) iff any minimizer T of has off- 

diagonal entry a ^ 0. 

It was shown in [22] that the question whether the system is superconducting can 
be reduced to the following linear criterion, which we will use heavily. (In [22], the 
results are proved for D = 3 and without the restriction to the reflection-symmetric 
subspace of L^(M^), but it was already observed in [16] that the statement holds as 
stated here.) 


Proposition 2.6 (Theorems 1 and 2 in [22]). Define the operator 

-pi 


Kx := 


tanh ( ^ 


as a Fourier multiplier and consider Kt + V in the Hilbert space 


L 


hD\ 


):={/£ LflR^) : /(x) = /(-x) o.e.}. 


sym \ 


( 2 . 6 ) 


(2.7) 


Then: 

(i) the system is superconducting in the sense of Definition 2.5 iff Kt + V has at 
least one negative eigenvalue. 

(ii) there exists a unique critical temperature 0 < Tc < oo such that 

Ktc + > 0 , 

inf spec(ibr + b^) < 0, VT < Tc. 


( 2 . 8 ) 


Tc is unique because the quadratic form associated with Kt is strictly monotone in 
T. In a nutshell, the reason why the operator Kt -|- V appears, is that it is the Hessian 
of the map 

° 0 )) 

at f = 0 with To the normal state of the system, see (2.12), and naturally, the positivity 
of the Hessian is related to minimality. For the details, we refer to [22]. In the following, 
we make 

Assumption 2.7. V is such that K > 0. 

By Theorem 3 in [22], H < 0 and V ^ 0 implies Tc > 0 in ZI = 3 and this result is 
stable under addition of a small positive part. 

Definition 2.8 (Ground-state degeneracy). We set 

n := dimker(Tr7’^-|-P). (2.9) 
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Remark 2.9. (i) We always have n < oo. The reason is that, by Assumption 2.2 on 

V, the essential spectrum of Kt + F is contained in [2r, oo). Therefore, zero is 
an isolated eigenvalue of Kt^ + F of finite multiplicity and so n < oo. 

(ii) A sufficient condition for n = 1 is that F < 0 and V ^0 [15, 24]. 

(hi) For Schrodinger operators —+ V, the ground state is non-degenerate by the 
Perron-Frobenius theorem. That is, one always has the analogue of n = 1 in that 
case. One may therefore wonder if n > 1 ever holds. In part III, we present a 
class of radial potentials such that for open intervals of parameter values, we have 
n > 1. In fact, one can tune the parameters such that ker(iPj'^ -|- V) lies in an 
arbitrary angular momentum sector. 

2.2 GL theory 

In GL theory, one aims to find “order parameters” that minimize the GL energy. The 
minimizers then describe the macroscopic relative density of superconducting charge 
carriers, up to spontaneous symmetry breaking. Microscopically, they describe the 
center of mass coordinate of the Cooper pair wave function a. In our case, translation- 
invariance implies that the order parameters are complex-valued constants, which are 
non-zero iff the system is superconducting. 

When n = 1 (and the system is translation-invariant), there is a single order pa¬ 
rameter '0 G C and for T < Tc the GL energy is of the all-familiar “mexican hat” 
shape 

T'^'^(0) = c|0|^ — djV'P, c, d>0. (2-10) 

Below, in Theorem 2.10, we show that for n > 1, the GL energy is of the form 

£'^^{a) = j /4(p)|a(p)|'^dp-y /2(p)|a(p)pdp (2.11) 

and 0 varies over the n-dimensional set ker(iL 7 ’^ -|- V). The functions and /2 are 
explicit; they are radial {p = jpj) and positive for T < Tc- 

Thus, we see that the mexican hat shape is characteristic for the translation- 
invariant case, even in the presence of degeneracies. However, there exists nontrivial 
coupling (i.e. mixed terms) between the different basis elements of ker(iL'r^ -|- V) in 
general. 


2.3 Result 


We write Fq for the minimizer of the free energy as in (2.2) but with V = 0. 

That is, Fq describes a free Fermi gas at temperature T and for this reason we call 
Fq the “normal state” of the system. From the Euler-Lagrange equation, one easily 
obtains 


fo(p) 


70 (p) 0 \ 

0 1 -7o(p) ) ’ 


( 2 . 12 ) 


where 

1-hexp((p2-^)/T) 


(2.13) 
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is the well-known Fermi-Dirac distribution. (Of course, Fq depends on ^ and T, but 
for the following we implicitly assume that it has the same values of //, T as the free 
energy under consideration.) 

We now state our first main result. It says that an appropriate n-component GL 
theory arises from BCS theory on the macroscopic scale and for temperatures close to 
Tc. Recall that p = |p|. 


Theorem 2.10. Let V satisfy Assumptions 2.2 and 2.1 and let p G R, T < Tc. 
Recall that n = dimker(R'r^ -|- V). Then: 

(i) AsTf Tc, 


2 


\ Tc J a&.er{KTc+V) 


(2.14) 


where is defined by 

yCL/ ^ [ 9i{{p^ - lj)/Tc) 

Tc Jrd (p2 - ^)/Tc 

1 r 1 


G JRO 


\KtM\^ |ci(p)l^dp 

l-fi^T,(p)Ha(p)Pdp. 


(2.15) 


2Tc 


Here we used the auxiliary functions 

tanh( 2 ;/ 2 ) 

9o{z) : = -^- 

9i{z) ■■= - 9oi^) = z~^9o{z) 


Kt{p) := 


p^ - p 


tanh ( 


1^-1 1 
2 cosh^(z/2) 


(2.16) 


(a) Moreover, ifT is an approximate minimizer of in the sense that 

J-#°yr)-J-#^"(r„)=(F^)V min f“(n) + e), (2.17) 

\ Tc J \a&.er{KTc+V) J 

for some 0 < e < M, then we can decompose its off-diagonal element a as 

a(p) = ^ ao(p) -h (2.18) 

where ||.^||2 = Om (2c — T) and Uq G ker(R'rc + k^) is an approximate minimizer 
of the GL energy, i.e. 

£^^{ 00 ) < min<5'^^ e -h Om {F - T). 

(Here Om means that the implicit constant depends on M.) 


11 



The idea is that near Tc, where superconductivity is weak, the normal state Tq is the 
prime competitor for the development of a small off-diagonal component a of the BCS 
minimizer. Theorem 2.10 then says that the lowest-order deviation from the normal 
state is well-described by a GLn whose coefficients are given explicitly as integrals over 
microscopic quantities. 

Remark 2.11. (i) We can equivalently rewrite the GL energy in terms of “order pa¬ 

rameters” V’l) • • •) V’n as follows. We fix an orthonormal basis {aj} of ker{KT^ + V) 
and decompose o G keic{KT^ + V) as a(p) = The basis coefficients 

■01, • •., V'n £ C are the n order parameters, each one corresponds to a different 
“type” of superconductivity Uj. The GL energy (2.15) can then be rewritten in 
the equivalent form 


T (V’l) • • • ) 0n) — ^ ^ dijipi'lpj. 




Here the “GL coefficients” Cijkm->dij are given by 


^ijkm — 


J_ f 9i{{P^ - r)/Tc) 
Tc Jro p‘^ - R 


\KtAp)\ aj(p)%(p)afc(p)am(p)dp 


~ 2Tr 


cosh2 ( ^ 


I^t(p)| ai(p)%(p)dp. 


(2.19) 


( 2 . 20 ) 

( 2 . 21 ) 


The minimum in (2.14) turns into the minimum over all 01,..., 0^ G C. In part 
n, we compute the integrals (2.20),(2.21) for special symmetry types and study 
the resulting minimization problem given by (2.19). 

(ii) If one assumes n = 1, this result is a corollary of Theorem 1 in [16], which is 
obtained by restricting it to translation-invariant systems. (When comparing, 
note that [16] rescale the BCS free energy to macroscopic units.) In this case, the 
microscopically derived GL theory is simply of the form (2.10). 

(hi) Note that the error term in (2.14) is 0{Tc — T) higher than the order at which 
the GL energy enters. Such an error bound is probably optimal because the 
semiclassical expansion of Lemma 5.4 will contribute terms at this order. It 
improves on the error term that one would obtain from Theorem 1 of [16] in the 
case n = 1. 

We note that writing minT^^ in the above theorem is justified because 

Proposition 2.12. The microscopically derived Ginzburg-Landau energy satisfies infc" > 
—oo. Moreover, the infimum is attained. 

When T > Tc, it was proved in [22] that the unique minimizer of T i—)■ is 

the normal state Tq. In other words, the left-hand side in (2.14) vanishes identically 
for all T > Tc- Nonetheless, one can still ask if GL theory describes approximate 
minimizers of the BCS free energy similarly to Theorem 2.10 (ii) when T — Tc is 
positive but small. Indeed, above Tc approximate minimizers must have small GL 
order parameters (as one would expect): 
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Proposition 2.13. Suppose T > Tc and T satisfies 

with 0 < £ < M. Let {aj} be any choice of basis for ker{KT^ + V). 
Then, there exist V'l, • • •, V’n £ C” and ^ G such that 


Ti 2 — a = 



with ||.f II 2 = Om {T - Tc) and 

n 

Y.\A\^ <^ + OMiTc-T) ( 2 . 22 ) 

1 ^min 

2=1 

as T Tc- Here Amin > 0 is a system-dependent parameter. 

3 Part II: Examples with d-wave order param¬ 
eters 

3.1 Angular momentum sectors 

In order to explicitly compute the GL coefficients given by formulae (2.20), (2.21), we 
make some assumptions on the potential V. First and foremost, we assume that V is 
radially symmetric. We can then decompose L^(M^) into angular momentum sectors. 
We review here some basic facts about these and establish notation. For the spherical 
harmonics, we use the definition 


F)-(^,^) = (S + m)l Pr{cos^)e^-^, (3.1) 

where P™ is the associated Legendre function, which we define with a factor of (—1)™ 
relative to the Legendre polynomial Pm- While we will use the in the proofs, it 
will be convenient to state the results in the basis of real-valued spherical harmonics 
defined by 



^ (y- - (-l)-yr-), 

if m < 0 



yO 

-'0 ) 

if m = 0 

(3.2) 


i {vr + i-iry,-"') . 

if m > 0. 


span{Y[^}m=- 

— span{y(^m}m=— 

and define 



= L^(M+; r^dr) (g) Si, 

(r = |x|). 

(3.3) 
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We employ the usual physics terminology 


"Ho = {s-waves}, T-ii = {p-waves}, 'H 2 = {d-waves}. (3-4) 

2 2 

Note that %q is just the set of spherically symmetric functions and > 2,2 oc is the 

dx 2 _y 2 -wave in this classification. In analogy to Fourier series, we have the orthogonal 
decomposition [45] 

00 

l2(m3) = 0^,. (3.5) 

1=0 

Recall that r = |x|. The Laplacian in 3-dimensional polar coordinates reads 


V2 

y/S _ y72 I 


(3.6) 


where = r ‘^dr{r‘^dr) and V ^2 is the Laplace-Beltrami operator, which acts on 
spherical harmonics by 

-Vl2Yi,m = lil + m,m. (3.7) 

Since Kt commutes with the Laplacian and V clearly leaves the decomposition (3.5) 
invariant, we observe that the eigenstates of Kt + V can be labeled by I (in physics 
terminology, I is a “good quantum number”). To make contact with unconventional 
superconductors, we will suppose we are in either of the two cases: 

• kei{KT^ + R) = span{/92} ® S 2 , “pure d-wave case” 

• ker(iL'r^ + R) = spanjpo} ® 5o + span{/92} ® ^ 2 , “mixed (s -j- d)-wave case”. 

Here po,p 2 £ L^(]R_i_; r^dr) are radial functions. They are determined as the ground 
states of an appropriate /-dependent operator acting on radial functions. We assume 
that these radial ground states are non-degenerate for simplicity. This assumption is 
satished for the examples we give in part III, but may not be satisfied in general. 


3.2 Results 

3.2.1 The pure cZ-wave case in three dimensions 

Theorem 3.1 (Pure d-wave case, 3D). Let D = 3. Let V be such that Theorem 2.10 
applies and such that ker(iL'r^ + V) = span{/92} ® ^2 for some 0 ^ p 2 £ L^(M+; r^dr). 
Let {a2,m}m=-2,...,2 be an orthonormal basis of the kernel such that 

a 2 ,m{p) = Q{p)y2,7n{'&,T) (3-8) 

for an appropriate g £ L^(M+;p^dp) (explicitly, g is the Fourier-Bessel transform (7.6) 
of g). Let tprn denote the GL order parameter corresponding to 0 , 2 ,m for —2 < m < 2. 
Then: 


(i) The GL energy that arises from BGS theory as described in Theorem 2.10 reads 


cGL 
^d-wave 


(V’-2, ...,'lp2) = 


5c 

Ilvr 


j] IV'n 


- T -T + 


\m =—2 


m =—2 


. (3.9) 
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where r := ^ and 


c= h{p)dp, d= f 2 {p)dp. 
Jo Jo 


Here, we introduced the positive and radially symmetric functions 


(3.10) 


U{p) = 


p 


,2 51 ( 4 ^ 


Tf p^-p 


h{p) =: 


' c cosh^ 


\KTAp)Q{pt 


\KtSp)q{p)? 


2Tc 


See (2.16) for the definition of gi and Kt{p)- 
(a) We have ~ “ rninimizers is 


(3.11) 


Md-wave = I (V’-2,..., V’ 2 ) G ^ IV’ml^ = T and ^ V’m = 0 4 (3-12) 


m =—2 


m =—2 


Remark 3.2. (i) The existence of V such that the assumption on ker(i^'r^ + V) holds 

for an open interval of parameter values follows from statement (i) of Theorem 
4.1 by choosing Iq = 2. 

(ii) Observe that the minimization problem in (3.9) is trivial, i.e. (ii) is immediate. 

(iii) Recall that we normalized the GL order parameters such that they are related 
to the Cooper pair wave function via (2.18). For the special case (3.12), we 
see that a minimizing vector will have absolute value ^/T. We can then reduce 
to the case where a minimizing vector lies on the unit sphere by rescaling the 
order parameters. The advantage of this other normalization is that it allows 
to interpret the absolute value of the order parameters as relative densities of 
superconducting charge carriers. 

We discuss what symmetry of one can expect. First of all, GL theory always 
has the global U{1) gauge symmetry ifj i—)■ (this is due to the presence of the 

absolute value signs in (2.19)). Second, 50(3) acts on spherical harmonics by pre¬ 
composition, i.e. for g G 50(3) and u) G 

m' 

where G 0(21 + 1) is the analogue of the well-known Wigner d-matrix for real 
spherical harmonics [3]. By changing the angular integration variable in (2.20) and 
(2.21) from giv to uj, it is easy to see that 

where we introduced if = {ip- 2 , ■ ■ ■ ,'02)- Since is reflection-symmetric for even I, 
we can extend the action to all of 0(3) and retain the invariance of This shows 
that we can expect to have symmetry groups 0(1) and 0(3). However; 
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Corollary 3.3. In the situation of Theorem 3.1: 

(i) For all f G [0, 27r), IZ G 0(5) and G C®, 

^ (3.13) 

Moreover, 0(5) acts transitively and faithfully on Aid-wave- 
(a) Aid-wave is a 7-dimensional manifold in 

(in) Any minimizer of has at least two non-zero entries ipj. 

Remark 3.4. (i) Surprisingly, the emergent symmetry group 0(5) is considerably 

larger than the 0(3)-symmetry discussed above. (Recall also that from above 
is in 0(5), so that the 0(3)-symmetry is really contained in the 0(5)-symmetry.) 
The particularly nice form of the 0(5) action is a consequence of choosing the 
real-valued spherical harmonics as a basis. 

(ii) We interpret faithfulness of the group action as saying that Af^-wave is “truly” 
invariant under the full 0(5). 

(hi) Transitivity means that the set of minimizers Aid-wa.ve is a single orbit under 
the 0(5) symmetry. In other words, there exists a unique minimizer modulo 
symmetry. 

(iv) We interpret (hi) as a proof of non-trivial coupling between the real-valued d-wave 
channels (it is of course a basis-dependent statement). 

Proof. The invariance under multiplication by e*'^ is trivial. To see the 0(5) symmetry, 
we use real coordinates because they also provide an interesting change in perspective. 
Writing if = x iy with x, y G M®, the GL energy becomes 

£^^{x + iy) = ^ (^(x^ -l-if - rf - ^ + \x ■ . (3.14) 

This is clearly invariant under the 0(5)-action x -\- iy TZx -\- iIZy. We can rewrite 

the set of minimizers as 

Aid-wa.ve = |(x,y) G X : x'^ = jf = X ■ y = o'^ . (3.15) 

Without loss of generality, we may set r/2 = 1, so that Aid-wiwe is just the set of pairs 
of orthonormal M^-vectors. To see that the 0(5)-action is transitive, consider the orbit 
of (61,62) G Aid-wuve, namely {(TZei,TZe 2 ) ■ TZ G 0(5)}. Since any two orthonormal 
vectors can appear as the first two columns of an orthogonal matrix, we have transi¬ 
tivity. To see that the action is faithful, note that for any two distinct IZ,IZ £ 0(5), 
there exists e* such that TZci / IZei. 

For (ii), we employ the implicit function theorem and observe that the Jacobian as¬ 
sociated with the functions x^,if‘,x-y from (3.15) has rank 3. Finally, (iii) is immediate 
from (3.12). □ 
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3.2.2 The pure d-wacve case in two dimensions 

Note that the two-dimensional analogue of the space Si, namely the homogeneous 
polynomials of order / on is spanned by cos{lip) and sin(/(/9). Thus assumption (3.16) 
below is the two-dimensional analogue of the assumption kei{KT^ + V) = span{/92}®‘52 
in Theorem 2.10 above. 


Theorem 3.5 (Pure d-wave case, 2D). Let D = 2. Let V be such that Theorem 2.10 
applies and such that keic{KT^ + V) = sY>aii.{axy,a^ 2 _y 2 } with 




(P) = 0(P) 


cos(2(/9) 


TT 


dxy{p) = Q{p) 


sin(2(/9) 


TT 


(3.16) 


for an appropriate, normalized 0 ^ g ^ L^(M+,pdp). Let ’’Pxy denote the 

corresponding GL order parameters. Then: 

(i) The GL energy that arises from BGS theory as described in Theorem 2.10 reads 

^d-wave,2Di'^x^—y‘^} f^xy) 


C 





1 

+ 2 


1^12 _y2 + 


xy 


(3.17) 


where c, d are defined in the same way as in Theorem 3.1 with f 2 {p), fiip) replaced 
by f2ip)/p,hip)/p- 

(a) We have aiin 2 _d = '^b.e set of minimizers is 


Md -wave^2D 

= ^{iPx^-yLi^xy) G : \'lpx^-yA‘^ + I'lpxyl^ = V’xy = o| 

= |(V',±#) G : IV'P = 


(3.18) 


Remark 3.6. (i) Statement (i) directly implies the first equality in (3.18) and the 

second equality is elementary. Note that the result can be conveniently stated in 
terms of the complex-valued spherical harmonics as well. 

(ii) From the second equation in (3.18), we see that the minimizers of the GL2 for a 
pure d-wave superconductor in two dimensions (in the cosine, sine basis) form a 
double cover of the minimizers of the usual “mexican-hat” GLl. 

(hi) A similar result holds for any pure angular momentum sector in two dimensions. 


3.2.3 The mixed (s -|- d)-wave case 

We write for the real part of a complex number z. 

Theorem 3.7 (Mixed (s-l-d)-wave case, 3D). Let D = 3. Let V be such that Theorem 
2.10 applies and such that ker(Arj'^ -|- F) = spanjpo} ® 5o -|- span{p 2 } ® S 2 for some 
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0 ^ Pq-,P 2 G r^dr). As an orthonormal basis, take 02 ,m as in Theorem 3.1 and 

as with 

as{p) = Qsip) t)- (3.19) 

Let 'ipm^ = —2,..., 2) and V's denote the GL order parameters corresponding to the 
respective basis functions. Then: 


(i) The microscopically derived GL energy reads 

= £?-Lve{'4^s) + £d-wave{i’-2, • • • , ^^2) + Scouplingi'^s, V’-2, • • • , ^^2) 

where T^-wavei'^-^^ • • •, "02) is given by (3.9), 


(3.20) 




.(4s) 

dvr 


-Ts) - 


T. 


(3.21) 


with Ts = ^ j and 


^couplingi^s 1 ' 0 — 2 ; . . . ; ' 02 ) 


.(2s) 


27r 

\ 


2|0.P ^ + K 


m =—2 


0s^ ( ^ V’( 


\m=—2 


+ 


Ttt 




0s(20o|0Op+ X] l"^l(-l)™’^H20o|0mP + 000^ 

m=itl,2 


+ \/3K 


0s ^ m (20210m P + 020m) 


m=±l 


+2\/3K + 2V^_2K [V’l^-i])] j . 

(3.22) 

T/ie coefficients c,d are given by (3.10). Moreover, form = 1,2,4, we introduced 


roo roo 

C^ms)= f4{p)gs{p)'^ dp, (i(2s) ^ / f2{p)gs{pf dp, 
Jo Jo 


with 02, 04 as in (3.11) and 

9s{p) 


Qs{p) 

q{p) 


(3.23) 


(3.24) 


(ii) The following are equivalent: 

• for all sufficiently small s > 0, and for any minimizer (0_2,...,02) of 
^d-wave’ Ihcrc cxists ips with 10s I = e such that 


^5 + d)-wavei'^s,i^-2, . . . , 02 ) < Sitavei'^-2, • • • , 02 ) 


(3-25) 
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(Hi) The following are equivalent: 

• for all sufficiently small e > 0, and for any minimizer ips of Sf.lfaveJ there 
exists {ip- 2 , • • •, ^ 2 ) with |'0m| < e for m = —2,..., 2 such that 

^5 + d)-wavei'^s, f-2, • • • , ^^ 2 ) < S?-Lve(^s) = (3.26) 

We see that + £;)_wave yields a much richer GL theory than Especially 

the terms which depend on the relative phases of several GL order parameters make 
this a rather challenging minimization problem. Accordingly, we no longer have an 
explicit characterization of the set of minimizers. However, using (ii) and (iii) above, 
we immediately obtain 

Corollary 3.8 (Non-trivial coupling of s- and d-waves). In the situation of Theo¬ 
rem 3.7 suppose that and < dc^^®\ Then any minimizer 

(fs, f-2, ■■■,^ 2 ) of £^/f d)-wave Satisfy and V'm / 0 for some -2 <m <2. 

Remark 3.9. (i) The existence of V such that the assumption on ker(iLj'^ + V) in 

Theorem 3.7 holds for appropriate parameter values follows from statement (ii) 
of Theorem 4.1. 

(ii) Using the same method and the two-dimensional analogues of all quantities above, 
one can also compute the GL3 that arises for a two-dimensional isotropic (s -|- d)- 
wave superconductor 

+ d)-wSLve,2Di''Ps,fx^-y ^) i^xy) 

= Sc\Tp^2_y2\‘^ + Sclfxyl"^ + -T 2dR[f^2_y2‘^fly] 

+ 4.c\f^2_y2\‘^\fxy\‘^ + 4c(2^^3^[V’/(V'x2-y2 + i’ly)] 

-h 8c(^®)|V’sP(|V'a;2-y2p -h Ifxyl"^) 

- And {\fx 2 -yA^ + Ifxyl^) - dyrd^^^^lV’sP 

Its complexity lies somewhere between the GL theories in Theorems 3.5 and 3.7. 
Setting fxy = 0 (that is, we forbid the dxy channel ad hoc), we obtain the GL2 

+ d)-wa.\e,2D^'^Si i 0) 

= 3c\fx2-y-2\* + -h 4c(^^)K[V's^V’x2-y2] (3-28) 

82c*'^®^|V'sP|V'a;2-y2p — 47rd|V'j,2_j^2p — 47rd*'^^^ I V's P 

Compare this with T^wave 2 D Theorem 3.5. While one cannot complete the 
square because the coefficients differ in a way that depends on the microscopic 
details, notice that the only phase-dependent term is of the form 

4c(2®)3fi;[^V2^_^2] (3.29) 

with >0. It is then clear that for minimizers, the dx 2 _y 2 - and s-wave order 
parameters must have a relative phase of ±L 
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4 Part III: Radial potentials with ground states 
of arbitrary angular momentum 


In this part, D 


3 and ;U > 0. Recall that 


Kt{p) 


- P 

tanh 


(4.1) 


and the operator Kt is multiplication by the function Kt{p) in Fourier space. Recall 
the definition (2.5) of the Dirac delta potentials 


^a,r(x) = -A(I(|x| - R), 

for A, R > 0. 

The following theorem says that, given a non-negative integer Iq, we can choose 
parameter values for A, R from appropriate open intervals such that the zero-energy 
ground state sector of Kt^ + V\r lies entirely within the angular momentum sector 

Theorem 4.1. (i) Let Iq be a non-negative integer. For every R > 0, there exist an 

open interval I C M+ and A* > 0 such that for all p G I and all A G (0, A*) there 
exists Tc > 0 such that 


inf spec{KT, Vx,r) = 0, 

ker{KT, + Vx,r) = spanjpij (g) Si^, 
inf spec(Rrr -|- R) < 0, VT < Tc. 


Explicitly, the (non-normalized) radial part is 

Jl,+ drp)\+i{Rp) 


Ploir) = r 


- ^- 1/2 


P- 


KtSp) 


dp. 


(4.2) 

(4.3) 

(4.4) 


(4.5) 


(ii) For every R > 0, there exists T* > 0 such that for all Tc < T^, there exist p,X > 0 
such that 


inf spec{KTc + Vx,r) = 0, (4.6) 

ker{KTc + Vx,r) = spanjpo} ® 5o -G span{/92} ® S2, (4.7) 

inf spec(R:r + R) < 0, VT < Tc. (4.8) 

with po,P 2 as in (4.5). 

Remark 4.2. (i) To be completely precise, in (i) there exists Tq such that the analogue 

of (4.2)-(4.4) holds with Tq in place of Tc. Then Tc = Tq by definition (2.8). 

(ii) The parameter R can be removed by rescaling p, A and T appropriately. 

(hi) In statement (i), for given ;U G /, A G (0, A*) and R > 0, Tc is given as the unique 
solution to the implicit relation 
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(iv) The fact that statement (i) holds for open intervals of ^ and A values is to be 
interpreted as saying that the occurrence of degenerate ground states for Kt^ + 
^A, R is “not un-generic”. This may be surprising at first sight, because if one 
replaces Kt + by the Schrodinger operator —-|- V, the Perron-Frobenius 
Theorem (see e.g. [35]) implies that the ground state is always simple. 

(v) The proof critically uses that Kt{p) is small (for small enough T) on the set 

{p : = /u}. Note that this set would be empty for // < 0. 

(vi) It is interesting to compare Theorem 4.1 with Theorem 2.2 from [15] which char¬ 
acterizes the critical temperature in the weak-coupling limit A —)■ 0 through an 
effective Hilbert-Schmidt operator acting only on of the Fermi sphere. For 
radial potentials, [15] shows that ker(iF'r^ -|- P) C T-Liq for all sufficiently small A 
iff Iq is the unique minimizer of 

where ji{z) = \/^Ji{z) is the spherical Bessel function of the first kind. While 
our proof here will be independent of [15], one can take V = Vx^r in (4.10) to see 
that the key fact needed to prove ker(iF'r^ -|- P) C T-Iiq is that there is a point at 
which jf^ > supi^i^jf. This is the content of Theorem A.l. 

We conclude by discussing the conceivable extensions of Theorem 4.1. Statement 
(i) also holds if Kr -|- P is defined on all of L^(M^) instead of just on , so 

there is nothing special about even functions in (i). 

Statement (ii) can not be generalized as much: (a) it will not hold when odd 
functions are also considered and (b) it does not generalize to arbitrary pairs {Iq, /o + 2) 
with lo even. The reason is that, as demonstrated within the proof of Theorem 4.1, 
for small enough T, (ii) is equivalent to the existence of a point where ^ 1/2 > >Ji+if 2 
for all even / > 1. The generalizations to more /-values described above require the 
analogous inequalities for Bessel functions. However, these facts will not hold in the 
cases above, as becomes plausible when considering Figure 1. 


5 Proofs for part I 

The strategy of the proof follows [16]. 

We introduce the family of BCS states Fa from which the trial state generating the 
upper bound will be chosen. The relative entropy identity (5.6) rewrites the difference 
of BCS free energies as terms involving Fa. 

The main simplification of our proof as compared to [16] is then in the “semiclas- 
sical” Theorem 5.3. While [16] requires elaborate semiclassical analysis for analogous 
results, the proof in our technically simpler translation-invariant case reduces to an 
ordinary Taylor expansion. 

Afterwards, we discuss how one concludes Theorem 2.10 by separately proving an 
upper and a lower bound. In the lower bound, the degeneracy requires modifying the 
arguments from [16] slightly. 
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5.1 Relative entropy identity 

All integrals are over unless specified otherwise. We introduce the family of oper¬ 
ators 

A 
-f) 

Here A is an even function on and we have introduced 

(5.2) 

the energy of a single unpaired electron of momentum p. Note that the choice A = 0 
in (5.1) indeed yields the normal state Fq dehned in (2.12). 

Recall that F is a BCS state iff 0 < F <1 and F is of the form (2.1). 

Proposition 5.1. Fa defined by (5.1) is a BCS state and 



Fa := 


1 


1 -I- exp(iFA/7') 



rA(p) 


7a(p) aA(p) \ 
Sa(p) i-7A(p)y 


with 


7a(p) 

Sa(p) 

Ea{p) 


A(p) tanh(.gA(p)/(2^)) 
2 E^ip) 

VHp)^ + l^(p)P- 


(5.3) 

(5.4) 

(5.5) 


Proof. It is obvious from (5.1) that 0 < Fa(p) < 1. Since = E‘^12 and since 

tanh(a;)/x only depends on it follows that 


Fa =- ^ - = \(l- tanh(iFA/(2r 

1 + exp(iFA/T) 2 V 

^ ^ ~ ^tanh(F;A/(2r))^ , 


which yields (5.3) and (5.4). 


□ 


We now give an identity which rewrites the difference in terms 

of more manageable quantities involving Fa, one of them is the relative entropy. 


Proposition 5.2 (Relative Entropy Identity, [16]). Let F he an admissible BCS state 
and 0 G Set A = 2Eo. It holds that 


CS(F) _ j-"^’^'(ro) 

= -|/tVc. 

T 


log (l + - log (l + 


dp 


(5.6) 



where i3(r,rA) is the relative entropy defined by 

S){T,Ta) := I Trc2 [</>(?, f a)] dp. (5.7) 

Here we introduced 

(j){x, y) = x(log(x) - log(y)) + (1 - x)(log(l - x) - log(l - y)), VO < x, y < 1. 

Proof. This is a computation, see [16] or [17]. □ 

For the sake of comparability with [16], we note that in the translation-invariant 
case the L^-trace per unit volume of a locally trace-class operator (which they denote 
by Tr) is just the integral of its Fourier transform and so 

Tr[r] = / tVc 2 [r(p)l dp. 

Jro l j 

5.2 “Semiclassical” expansion 

We prove Theorem 5.3 by a Taylor expansion, which is sufficient because of the simpli¬ 
fications introduced by the translation-invariance. The analogous results in [16] require 
many more pages of challenging semiclassical analysis. 


5.2.1 The result and the key lemma 

Recall the definition of gi in (2.16). The following is the main consequence of the 
Taylor expansion 

Theorem 5.3. Let A = 2hVa for some o G ker(iVT'c + k). Define h > 0 by 



Then, as /i —)• 0, 

T-#^^(rA) - ^^(Fo) = h^E2 + 0(/r®), (5.9) 

where 


giim/Tc) 


it(p)rdp- — 


cosh2(f}(p)/(2rc)) 


l^(p)Pdp (5.10) 


with t = 2Va. 


We emphasize that this is the place where the effective gap operator appears in the 
analysis. The choice o G kei{KT,, + V) ensures that there are no 0(/i^) terms in the 
expansion (5.9). 

The theorem follows from the key 
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Lemma 5.4. Let Fa he given by (5.1) with A(p) = ht{p) for a function t satisfying 


t G with I 


q = oo if = 1, 
4<g<oo ifZl = 2, 
(g = 6 if = 3. 


Then, as h ^ 0, 

(^) 


-y/Trc. 


log (l + - log (l + 

= h^Ei + h^E2 + 0{h^) 

where E 2 is defined by (5.10) and (see (2.16) for pq) 


dp 


(5.11) 


(5.12) 


Ei = - 


ATr 


9o{Hp)/T c)\t{p)\‘^ dp, 


(5.13) 


(ii) 

\\aA - 4>\\m = 0{h^) (5.14) 

with 4>{p) = -h^i^go{t){p)/Tc). 

This may be compared to Theorems 2 and 3 in [16]. 

To conclude Theorem 5.3 from the key lemma, we need a regularity result for the 
translation-invariant operator. 

Proposition 5.5. Let a G satisfy {Kt,, + V)a = 0. Then, a G L°°(M^). Let 

t := Va and (p) := (1 -|-p^)^/^. Then, {p)~^t G LfifEP) and t satisfies (5.11). 

Proof. Recall Assumption 2.2 on the potential V. When V G L^^ then the result 
follows from Proposition 2 in [16]. For the potentials Va,r in H = 3, the regularity 
properties can be read off directly from the explicit solution of the eigenvalue problem 
(iFj-^ -|- Vx^R)a = 0, see (7.8) in the proof of Lemma 7.1 for its Fourier representation. 
Indeed, since and the Bessel function of the first kind fTi+ 1/2 are smooth and 
bounded with 77;+i/2(0) = 0 and since E < 2T, we get a G L°°. Moreover, 

+( \ V /Q N ^'+1/2(4*77) 

HP) « Vl^rn I'd, T) -^- 

y/p 

and since 77;_|_i/2 also decays like for large p-values, the regularity properties of 

t follow. In D = 1,2, one can again solve the eigenvalue problem {Kt^ + Vx^jfija = 0 
explicitly and obtains the claimed regularity by similar considerations. The details are 
left to the reader. □ 

Proof of Theorem 5.3. First, note that t = 2Va has all the regularity properties needed 
to apply (i), thanks to Proposition 5.5. We invoke the relative entropy identity (5.6) 
and use Lemma 5.4 to find 

-F#^^(Fa) - ^^(Fo) 

= h^Ei + h^E2-h^ [ Pja 


Mx-|- J V\aA — ha\^dx + 0{h^). 


(5.15) 
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Observe that 


g^{i){p)/Tc)=T,K^l{p). (5.16) 

By Plancherel and the eigenvalue equation {Kt^ + V)a = 0, (5.15) becomes 

-^r^^(rA) - -F'F^^(ro) = h^E2 + j P|aA - ha\^dx + 0(h®). 

Thus, it remains to show 

J P(x)|aA(x) —/ia(x)p dx = 0(/i®). (5-17) 

To see this, recall that V is form-bounded with respect to —V^, so it suffices to prove 
that lloA — ha\\}ji = 0{h^). Using the eigenvalue equation and (5.16), 

a(p) = -K:^^{p)Va{Y>) = -^5o(f}(p)/7):) 
and so (5.17) follows from Lemma 5.4 (ii). □ 

5.2.2 Proof of Lemma 5.4 

Proof of (i) We have 

log (l + = -#A(p)/(2r) + logcosh(#A(p)/(2r)). 

Observe that Tr £2 77a (p) = 0, that x i—)■ coshx is an even function and that 

-f^A(p)^ = 77 a(p)^ 72. We find 

Trc2 log(l + e~^^/'^^ = Trc2 logcosh(#A(p)/(2r)) 

= 21ogcosh(£'A(p)/(2r)). 

This and a similar computation for A = 0 show that 

- ^ y Trc 2 log (^1 -7 - log (^1 -7 dp 

= —T J (logcosh(i7A/(2T)) — logcosh(l)/(2T))) dp. (5.18) 

We denote the function in (5.18) by 

/('*') (1^) - (2?^)))' 

where we wrote E{h?) for E/\ and T{h?) = Tc(l — h?). Note that E' = \t\^/{2E) and 
recall the definition (2.16) of qq and gi. By an easy computation 

/( 0 ) = 0 , 

/'(0) = -5o(f)/Te)^, 

1 .. . ^ gi{^/Tc) _ 1 

2-' ^ ^ f) I6r2 cosh2(f)/(2re))8rc' 
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With this, we can expand (5.18) as follows 


Tr, 


'£2 


loghA - log To I dp 
2 


(5.19) 


= i/o((i/rc)^dp 

1 /■<7i(f)/Tc),„4 




i6r2 






cosh^ I ^ 


'dp +0(/i®). (5.20) 


It remains to check that the 0{h^) term is indeed finite. Using the Lagrange remainder 
in Taylor’s formula, it suffices to show 


sup ^|/'"(5)| dp < oo. 
0<<5</i2 


(5.21) 


We will control this quantity in terms of appropriate integrals over t which are finite 
by our assumptions on t. We introduce the function 


®(*) = icoJ(z/2) 


(5.22) 


By a straightforward computation 


1 


r{6) 


8T{6)^ [l2E{5y V E{5)/T{5) 


^g^{E{5)/T{5)) 

3 - g 2 [E{d)/T{d)) 


Tr 


2T{5) 


\t\%{EiS)/Ti6)) 


+T-. 


■2ui2 


COsh2 ^ 
cosn I 2r(5) 




Note that, for h? small enough, Tc/2 < T{5) < 2Tc for all 0 < 5 < h'^. Using this and 
the fact that and g 2 {z) are monotone decreasing for z > 0, we can estimate 


[ sup I/'"(5)I dp 

J 0<(5</i2 u! 


0<(5</i2 
<^1 


I® sup E{S)-‘^ 

0<5</i2 


8^-mm-g.m)/Tis)) 


m/m 


dp 


+ C 2 J\t\^g2ii)/m))dp 


+ C 3 


1 


cosh^(f)/(4T,)) 


+ £/2(f)/(2rc)) sup U((J)^ dp 


0<<5</i2 


(5.23) 

(5.24) 


Here Ci, 6*2, C3 denote constants which depend on D,Tc and may change from line to 
line in the following. For definiteness, assume D = 3. The arguments for 44 = 1, 2 are 
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similar. Since g 2 {z) is a bounded function that decays exponentially for large z, we 
can use Cauchy-Schwarz and the fact that i){p) ~ C{p)‘^ for large p to conclude 

C2 j |t|V(f)/(2Tc))dp < (72 j {\tf + dp 

and the right-hand side is finite by Proposition 5.5. Using that < 

for small enough h, the same argument applies to the C 3 term in (5.24). 

The Cl term in (5.24) contains a factor E{d)~‘^ which looks troubling because, as 
(5 —)■ 0, it is of the form and thus singular on the sphere {p : p^ = ;u} if ^ > 0. 
For the radial integration, this singularity would not be integrable (and we have not 
even considered the factor |t|® yet). However, the singularity is canceled by the factor 
3gi{z)jz — g 2 {z) with 2 ; = E{5)/T{S) in (5.24). To see this, recall the definition (2.16) 
of gi and (5.22) of g 2 and observe that gi{z)/z and g 2 {z) are both even functions. 
Using the power series representation for and tanh, it is elementary to check that 
in the expansion of 3gi{z)/z — g 2 {z) the coefficients of order z~‘^ and z^ vanish and so 
the lowest non-vanishing order is Therefore, the singularity is removed and since 
gi{z)/z and g 2 are bounded, we get 


sup E{6) ^ 

0<<5<h2 


^ gi{E{d)/T{5)) 

m/m 


g2{E{5)/T{5)) 


< C < 00 . 


Since f |t|® dp < 00 by our assumption on t, the (7i term in (5.24) is finite and we have 
proved (5.21). 

Proof of (ii) From (5.4) we have 


Sa(p) 


-h 


i(p) 

2 T 


go(EA(p)/T). 


Therefore 

||aA - <^11^1 = - /( 0 )P dp, 

where we introduced the function 

2s _ gojEm/m)) 

2T{h^) 


(5.25) 

(5.26) 


Recall that gQ = —gi- Using this and the fact that for h‘^ small enough, Tc/2 < T{5) < 
2Tc for all 0 < (5 < /i^, Taylor’s theorem with Lagrange remainder yields 

< Ch^ {\gii{E{ 6 )/T{ 6 ))\ + \gi{E{ 6 )/T{5))\ +U(5))) . 


Note that go{z) and g\{z)lz are monotone decreasing and so 

|'7i(f)/(2Te)) 


\f{h^) - /(0)| < Ch^ \go{i)/m))\ + C" sup 

0<<5<h2 




{\t\^ + EiSf) 


< 


Ch^ (\go{h/{2T,))\ 


+ 




d 


(|tP + d^) 


<ch^{\t\m+m) 
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where in the second step we used that E{5) = + |tp for small enough 

h and in the third step we used go{z) < C{z)~^ as well as gi{z)jz < C{z)~^. Assume 
D = 3 for definiteness. We can bound (5.25) as follows 

j{pf\t\Vih^)-f mi''dp <Ch^ I {\tf{p)~"^ + m{p)~^) dp = Ch\ 

where the last equality holds by the assumption on t. This proves (ii). □ 

5.3 Proof of Theorem 2.10 

We follow the strategy in [16]. That is, we prove theorem Theorem 2.10 (i) by separately 
proving an upper and a lower bound on the left-hand side in (2.14). The upper bound 
follows by choosing an appropriate trial state Ta and using the semiclassical expansion 
of the BCS free energy in the form of Theorem 5.3. For the lower bound, we show 
that the chosen trial states Ta indeed describe any approximate minimizer T to lowest 
order in h (this is precisely statement (ii) in Theorem 2.10) and conclude by using the 
semiclassical expansion once again. 

5.3.1 Upper bound 

Recall the definition of h in (5.8). In this section we prove 

minjf ^^(T) - U^^^(ro) < min £^^{a) + 0(/r®), (5.27) 

r oeker(ii'Tj.+V) 

where is given by (2.15). 

We get this by using the trial state Ta, defined by (5.1) with the choice 

A = 2h(^) (5.28) 

where o G ker(Ar 7 ’^ -|- V) minimizes (recall that minimizers exist by Proposition 
2.12). Then, (5.27) follows from Theorem 5.3 and the fact that evaluating the definition 
(5.10) of E 2 for the choice 

^(p) = A(p)//i = -2Kt,{p)^{p) 

produces the definition (2.15) of □ 

5.3.2 Lower bound: Part A 

Following [16], we will prove the lower bound in (2.14) in conjunction with statement 
(ii) about approximate minimizers. We consider any BCS state F satisfying 

jzBCS^Y) - 77^^^(Fo) < 0(/i^). (5.29) 

Note that we may restrict to such F when minimizing thanks to the upper 

bound (5.27) and that (5.29) still includes the approximate minimizers considered in 
(ii). In Part A, we prove Proposition 5.6, which says that the off-diagonal element ol 



of such a r will be close to a minimizer of . In Part B, we will use this to get 

> 0{h^) for A of the form (5.28) and hence 

jrBC5(r) _ > j-#^^(rA) - -^'F^^(ro) + 0 {h^). 

Since we know (V’l , • • •, '^n)+0{h^) from Theorem 5.3, 

this will imply both the lower bound in (2.14) and statement (ii) about approximate 
minimizers. 

In the remainder of this section, we will prove: 

Proposition 5.6. SupposeT satisfies (5.29) and let P denote the orthogonal projection 
onto ker(ii'p, + V) and let = 1 — P. Then, ||Pa ||2 = 0{h) and ||P-'“q ;||2 = 0{h?‘). 

This implies statement (ii) in Theorem 2.10 with uq = h~^Pa. The proof of Propo¬ 
sition 5.6 will use the following lemma, which bounds the relative entropy i^(r,rA) 
from below in terms of a weighted Hilbert-Schmidt norm. The result without the sec¬ 
ond “bonus” term on the right-hand side first appeared in [23], the improved version 
is due to [16]. 

Lemma 5.7 (Lemma 1 in [16]). For any 0 < T < 1 and T^^^ = (1 -|- exp(iL))“^, it 
holds that 


i5(r,r(^))> j Trc2 


H 


tanh(iL/2) 


(f _fW)S 


dp 


(5.30) 


+ 3 /Trc 2 


(f(i-r)-f(^)(i-f(^)))^ 


dp. 


Proof. By the identity (5.7) in [16] and Klein’s inequality for 2x2 matrices, (5.30) 
even holds pointwise in p. □ 


Here is a quick outline of the proof of Proposition 5.6: Following [16], we rewrite 
by invoking the relative entropy identity (5.6). Then, we bound 
the right hand side from below by {a, {Kt + K)a), which is therefore negative due to 
(5.29). Since Kt^ -|- K > 0 with a spectral gap above zero, this will allow us to conclude 
that the part of a lying outside of ker(Krj'^ -|- K) must be small, more precisely that 
jja —PQ;jj 2 = 0{hfi). To get that llPajj 2 itself is 0{h), we use the second “bonus” term 
on the right-hand side of Lemma 5.7. 


Proof of Proposition 5.6. Step 1: We first apply the relative entropy identity (5.6) with 
the choice a = 0 to get 


0{h^) > - T#^^(ro) = |io(r. To) + j Kjal^ dx. (5.31) 

Next, we use Lemma 5.7. To evaluate the resulting expression, note that 


-7^— 77 — .<2X2) 

tanh(77o/(2T)) 

f (1 - P) - fo(l - To) = ( 7(1 - 7) - %(1 - 70) - jap) I2x2, 
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are diagonal matrices. We obtain 


|i3(r,ro)> I (KT(-)(7-7o)'dp + i^r(-)l«l" 

2 \ 

+ ^ (7(1 - 7) - 70(1 - 70) - lap) j dp. 

We estimate the first term using Kt{p) > 2r and find the lower bound 

j + 2 T{j - ^ (7(1 - 7) - 70(1 - 70) - |ap)^^ dp. 

By 

(x(l - x) - y(l - y)f < (x - yf, VO < x, y < 1 
and the triangle inequality, we get the pointwise estimate 

^2(7 - ^ (7(1 - 7) - %(1 - %) - |ap)^^ > ^|ap. 

Going back to (5.31), we have shown that 


4T 

~5 



+ {a, {KT + V)a) <0{h^). 


(5.32) 


Step 2: Next, we replace Kt by Kt^ in (5.32) to make use of the spectral gap of 
Ktc + This is an easy version of what is Step 2 of Part A in [16], which is more 
involved because it also removes the dependence on the external fields A, W. For us, 
it suffices to observe that 


d 

dr 


Kt{p) 


J- _ ^(p)^ 

2r2 sinh^(f)(p)/(2r)) 


is uniformly bounded in p for all h small enough such that T > Tc/2. By the mean-value 
theorem, \\Kt — Kt^\\oo < 0{h?‘). Using this on (5.32), we find 


4r 

T 



+ (a, (ATtc + y)oi) < 0 {h^)\\a \\2 


+ 0{h^). 


(5.33) 


Let K > 0 denote the size of the spectral gap of Kt^ -t- V above energy zero. We write 
a = Pa + P-*-a. Using (LCtc + V)Pa = 0, we obtain 


— ||a||4 -f K||P‘''a||2 < 0(/i^)||a||2 + 0{h‘^). (5.34) 

t) 

For the moment we drop the first term on the left-hand side of (5.34) and use orthog¬ 
onality to get 

llP^alP < 0{h^){\\Pa\\l + llP^alli) + 0{h^) 

which yields 

||r^a|p < 0{h^)\\Pa\\l + 0{h^). (5.35) 

Thus, both claims will follow, once we show ||Pq ;||2 = 0(h). 
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Step 3: Here the degeneracy requires a slight modification. We now drop the second 
term on the left-hand side of (5.34) to get 

||S||4 < 0 {h^/‘^)\\a\\y‘^ + 0{h). (5.36) 

By orthogonality and (5.35), 

||S||4 < 0 {h^/‘^)\\Fa\\y^ + 0{h), (5.37) 

On the right-hand side of (5.36) however, the replacement of 3 by Pa requires more 
work. By the triangle inequality for || • ||4 and (5.35) 

||3||4 > II^IU - > ll^iu - ||P^||2^^||P^||^^ 

>||7W||4-0(3i/2)||^||f IiK.IIJ/ 2. 

We use P-^a = 3 — Pa and |3p < 7(1 — 7) < 1/4 pointwise to find ||P-*-a||oo < 
I -|- 11 Halloo- It is slightly more convenient to conclude the argument by choosing an 
orthonormal basis {aj} for ker(H'rc -|- V). This allows us to write 

n 

Pa = h (5.38) 

i=i 

By Proposition 5.5, ||a)||oo < C for all j and therefore ||Ha||oo < 0(3)1^^100- We have 
shown 

||3||4 > ||iw||4 - 0(31/2)||P^||1/2(1 + 

Combining this with (5.37), we obtain 

ll^lk < 0{h^/^)\\Fa\\y^{l + + 0{h). (5.39) 

It remains to bound ||Ha||4 from below in terms of ||Ha||2. Let i? > 0. We split the 
integration domain into {p < H} and {p > R}. Applying Holder’s inequality to the 
former yields 

||/W||i<OH^/2||^||2^;^2|^|^^ j |S^||^.|dp (5.40) 

{p>R} 

where O > 0 denotes a constant independent of h,R. Note that for all i,j, Cauchy- 
Schwarz implies |ai||aj| G L^(M^) and so for Hq > 0 large enough, 

j l«*ll%|dp<^- 

*3 {p>Ro} 

We recall (5.39) to find 

\\Fa\\l < 0(3)||iW||2(l + LIV-loo) + ^3V|L + 0(3')- 
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Since the {aj} in (5.38) are orthonormal, /i|V’|oo < ll-Pc^lb < \^h\'tp\oo- This implies 

+ ow) ml < 0{lMoo + 0(1). (5.41) 

Let h be small enough such that the 1/2 + 0{h) term exceeds 1/4. We conclude that 
iV'loo < 0(1). Since ||T’a||2 < v^^lV'loo, it follows that ||T’a||2 < 0{h) as claimed. □ 


5.3.3 Lower bound: Part B 


We use once more the relative entropy identity (5.6). Together with Lemma 5.4 (i) 
and the eigenvalue equation, we get 




BCS( 


(5.42) 


=m£^^{Pa) + l^{T,T/^) + j V\a- Pa\‘^dx + 0 {h^). 

We see that to prove the lower bound it remains to show 

|l3(r, Ta) + jv\a- Pa\^ dx = Ta) + j V\P^a\‘^ dx > 0(/i®). (5.43) 

By Lemma 5.7 and the fact that x i—)• x/tanh(x) is a monotone function that depends 
only on we have 


^f)(r,rA) f Trc2 


= -,jTrc. 
1 


(r-TA) 

(f-TA) 


Ha 


tanh \^Ha/ (2r 
Ea 


> 2 / Trc. 


tanh {Ea/ (2T)) 
(T - rA)iLT(r - f a)1 dp. 


(r-TA) 


(r-LA) 


dp 

dp 


Since Kt > 0, we have for every fixed (i.e. /i-independent) 0 < e < 1, 

1 


- / TVc 2 [(T - rA)iLr(r - Ta) 


dp 


> J iLr|S —CKApdp 

> [ iLr|Kpdp-2K [ Kr^ai Pa — a A ] dp 


> (1 — e) J KT\P-^a\^ dp — Ce J iLrlPa —OApdp 

> (1 - e) j Kt\P^\‘^ dp + 0 {h^). 

In the last step, we used Lemma 5.4 (ii) and Kx{p) < C{p)‘^ to get 


(5.44) 



Using these estimates on (5.43) and setting ^ := P-^a, we see that it remains to show 
that there exists an h-independent choice of 0 < e < 1 such that 

{ti{l-s)KT + V)0>O{h^). (5.46) 

Recall from step 2 of the proof of Proposition 5.6 that \\Kt — R^rdloo < 0{h?‘). Since 
also ||^||2 = 0{h?‘) by Proposition 5.6, we get 

(^, ((1 - e)KT + U)0 = (e, ((1 - e)KT^ + V)i) + 0 {h^). 

We claim that there exists a constant c > 0 such that 

{i,{KT^ + V)i)>c{i,KTA)- (5.47) 

Choosing e sufficiently small will then give (1 — s)Kt^ + U).^) > 0. Thus, it remains 
to prove (5.47). Since V- is infinitesimally form-bounded with respect to Kt^, we have 
for any 5 > 0 

(1 - 5)Kt, < Kt, -V-+Cs<Kt, + V + Cs (5.48) 

or 

Kt, < Cl {Kt, + V) + C2 . (5.49) 

Now, on ker(Rr 7 ’^ + V)-^, it also holds that Kt^ -|- U — k > 0 where k > 0 denotes the 
gap size. Thus, for all A > 0, 

Kt, < (Cl + \){Kt, + U) + C 2 - Ak, onker(i^T, + V)^ (5.50) 

and choosing A = C 2 IK, we see that (5.47) follows. This proves (i). 

Statement (ii) was proved along the way: Any approximate minimizer satisfies 
(5.29) and hence Proposition 5.6 implies that its off-diagonal part can be split into 
a = Pa + ^ with ||^|| = 0{h?). Since P is the projection onto ker(A' 7 .^ -t- U), Pa G 
ker(Ar 7 ’^-|-U). Moreover, ao = h~^Pa approximately minimizes the GL energy because 
the proof of the lower bound shows that for all P satisfying (5.29) (not just for actual 
minimizers), 

j-BCS(r) _ + 0(/i®). 

This finishes the proof of Theorem 2.10. □ 

5.4 Proofs of Propositions 2.3, 2.12 and 2.13 

Proof of Proposition 2.3. For the potentials, this is a standard argument combin¬ 
ing Holder’s inequality and Sobolev’s inequality. 

Consider the potentials (2.5), i.e. U(x) = —A(5(|x| — R) with X,R > 0. Let / G 
77^ (M^). We hrst consider the case D = 1. Then 

{f,Vf) = -\\f{R)\\ 

We apply the simplest Sobolev inequality 

/ CXD 

|u'(x)|dx, Vu G (5.51) 

-CX) 
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(which follows from the fundamental theorem of calculus) with the choice u{s) = /(s)^. 
By (5.51) and Cauchy-Schwarz, we get 

/ OO 1 

^|/(x)/'(x)|dx<e||/'||i + -||/||i 

for any e > 0. This proves the claimed infinitesimal form-boundedness of V when 
D = l. 

Let now D = 2,3. We have 

(/,!//) =-A [ R^-^\fiRw)\^da{uj), 

J§D-1 

where da is the usual surface measure on Observe that the inequality (5.51) 

implies 

poo 

2sup|ii(s)|< / |ii'(s)|ds, Vtt G 1 Lq^’^(M+). 

s>0 JO 

We use this with the choice u{s) = s^~^f{su})‘^, pointwise in to G and find 

[ R^-^\f{Rw)fda{u:) 

f rfD-i \ 

< ^ J y 2 ^\f{suj)f + s^ ^\f{su;)dsf{suj)\jdsda{uj). 

Consider the first term in the parentheses. We split the integration domain into s > 1 
and s < 1 and estimate s^~‘^ < in the first region. By applying Holder’s inequality 
in the second region, we get 

p poo p pi 

/ / s^~^\f{saj)fdsda{u)<\\f\\l+ / s^~‘^\f{su;)fdsda{u;) 

JS^-iJO J§D-1 Jo 

/ /.I x3/5 

<11/112+^ S^-^/^dsj ll/lli 
= \\ff2 + C\\f\\l 

where C is a finite constant, since D — 8/3 > —1. The norm is infinitesimally 
form-bounded with respect to —by the usual argument via Sobolev’s inequality. 

We come to the second term in (5.52) in parentheses. By Cauchy-Schwarz, for every 
e > 0, it is bounded by 

ef f s^-^\dsf{suj)\'^dsda{uj) +^WfWl- 
JsD-i Jq 4e 

The first term is the quadratic form corresponding to (the negative of) the radial part of 
the Laplacian, see (3.6). It differs from the full Laplacian by a multiple of the Laplace- 
Beltrami operator — V|o_i, i.e. a nonnegative operator. This implies infinitesimal 
form-boundedness when D = 2,3. □ 
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Proof of Proposition 2.12. Recall (2.15) 





9i{{p^ - p)lTc) 

r 1 


cosh^ 



l^r,(p)ha(p)|^dp 

r IKtSp)^ |a(p)Pdp, 


We denote the quartic term by ^(a) and the quadratic term by —R(o). Note that 
A,B > 0 whenever o is not identically zero. 

We use the basis representation of the GL energy mentioned in Remark 2.11 (i). 
That is, we fix a basis {oj} of keic{KT^ + V) and write a(p) = with 

(V'l,..., V’n) £ C". Then we write 

(•01,... ,'0n) = Tw, L > 0, w G 5(C’"), 


where S{£P) is the unit sphere in C”. It follows that 


inf £:^^(0 i,.. 


• ) V’n) 


inf inf 

a;eS(Cn)i>0 

inf inf (L'^Aioj) — L^Bioj)) 
a;eS(Cn)L>0^ ^ 

mf ^ , , 

a;eS(Cn) 4yl(a;) 


and since A, B are continuous functions which never vanish on the compact set S{C^), 
the last infimum is finite and attained. □ 


Proof of Proposition 2.13. The same argument that proves Theorem 2.10 (ii) applies 
for T > Tc and yields the same result with the sign of the |ap term in the GL energy 
(2.15) flipped. Consequently, the unique minimizer of the GL energy is a = 0. To 
see coercivity of the GL energy around this minimizer, we drop the quartic term and 
rewrite the the quadratic term as in the proof of Proposition 2.12 above. We get 

n 

T'^^(01, ...,'lpn)> |0iP 

i=i 


with 


Xmin ■= min — 
a;eS(C") 2Tc 


cosh^ 


2Tc 


\KtSpW 




ujjajip) 


dp. 


Note that Xmin 
compact set. 


> 0, since it is the minimum of a positive, continuous function over a 

□ 
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6 Proofs for part II 

6.1 Setting 

We use the formulation of GL theory from Remark 2.11(i). We compute the GL 
coefficients Cijkm and dij given by formulae (2.20) and (2.21). They determine the GL 
energy £ 2 Lve : C® ^ M via 

2 _ 2 _ 

{^-2, ■ ■ ■ ^ ^ 
i,j,k,m=—2 iJ=—2 

It remains to pick a convenient basis to compute (2.20) and (2.21). Since the Fourier 
transform maps Tii to itself in a bijective fashion, see e.g. [45], we can choose 

dmip) = Q{p)Y^{'d,ip), p = {p,p,(p), (6.1) 

for an appropriate radial function q. We will denote the GL order parameter corre¬ 
sponding to dm (in the sense of (1.1)) by 'ipm with — 2 < m < 2. (Note that we use 
the ordinary spherical harmonics (3.1) as a basis because it is more convenient to 
do computations, but our final result is phrased in terms the basis of real spherical 
harmonics (3.2).) 

With the choice (6.1), equations (2.20),(2.21) for the GL coefficients read 

= Ip-^h{p)Yi{d,ip)Y^{d,ip)Yi{d, ip)Yr{d, if) dp ( 6 . 2 ) 

dij = - jP~^f2{p)Y.^{'d, (p)Y^id, (p){p) dp, (6.3) 

where i,j,k,m = —2,..., 2 and we used the functions / 2,/4 defined in (3.11). Note 
that / 2,/4 are positive (since gi defined by (2.16) satishes > 0 ) and radially 

symmetric. 

6.2 Proof of Theorem 3.1 

While the radial integrals in ( 6 .2),( 6 .3) depend on the details of the microscopic poten¬ 
tial V through Q, the integration over the angular variables can be performed explicitly. 
Since the spherical harmonics form an orthonormal family with respect to surface mea¬ 
sure on we immediately get 

dij — dSij 

where d > 0 is the result of the radial integration in (6.3), i.e. 

poo 

d= f2ip)dp (6.4) 

Jo 

and this is the second relation claimed in (3.10). 
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i 

j 

k 

m 

^ijkm ' 287r 

2 

2 

2 

2 

lOc 

2 

1 

2 

1 

5c 

1 

1 

1 

1 

10c 

0 

2 

0 

2 

5c 

1 

1 

0 

2 

0 

0 

1 

0 

1 

5c 

-1 

2 

-1 

2 

5c 

0 

1 

-1 

2 

0 

0 

0 

0 

0 

15c 

1 

-1 

1 

-1 

10c 

2 

-2 

2 

-2 

10c 

0 

0 

2 

-2 

5c 

0 

0 

1 

-1 

-5c 

1 

-1 

2 

-2 

-5c 


Table 1: Non-trivial equivalence classes of Ginzburg-Landau coefficients in the pure d-wave 
case, c is defined as the result of the radial integration (6.5). Notice that the case i + j = 0 
behaves rather differently. This is due to the fact that the “pair permutation” and “pair sign- 
flip” symmetries fall together in this case. We keep the factor 5 to ensure better comparability 
with Table 2 later on. 


Next, we consider (6.2). Firstly, note that Cijkm is always proportional to the result 
of the radial integration in (6.2), i.e. 

roo 

c= U{p)dp (6.5) 

Jo 

and this is the first relation claimed in (3.10). 

It remains to compute the angular part of the integral in (6.2). We express the 
product of two spherical harmonics of angular momentum / = 2 as a linear combination 
of spherical harmonics of angular momentum ranging from / = 0 to / = 4. The general 
relation involves the well-tabulated Clebsch-Gordan coefficients, which we denote by 
{li,l 2 ]mi,m 2 \L; M), and can be found in textbooks on quantum mechanics (see e.g. 
[10] p. 1046): 


L=\h-l2\^ ^ ’ 


( 6 . 6 ) 


X {li,l2;mi,m2\L;mi -|-m2)T™^'''™'^(d, (/?). 


Physically, this corresponds to expressing a pair of particles, uncorrelated in the an¬ 
gular variable, in terms of a wave function for the composite system. Since the total 
angular momentum of the composite system is not determined uniquely by the product 
wavefunction on the left-hand side, the sum over L appears on the right. However, 
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the total z-component of the angular momentum is determined to be mi + m2. This 
“selection rule” will greatly restrict which Cijkm may be non-zero. 

Now, we can use the orthonormality of the spherical harmonics to compute the 
angular integrals and find 

_^ 25c 

Cijkm = 4 (2L + ll h + 3 ) 

L=o,2,4 ’ (6.7) 

X (2, 2, /t, m|Z/, k “t" TTl)Si-^j^k-\-m^ 

where we used that the Clebsch-Gordan coefficients are real-valued and that (Zi, I 2 ] 0,0|L; 0) 
0 unless L is even [10]. Note that the selection rule from above yielded the necessary 
relation i -|- j = A: -|- m for Cijkm / 0. 

There are further symmetries: Considering the original expression (6.2), we that 
Cijkm = Cjikm = Cijmk- Since (6.7) shows Cijkm £ 1^, (6.2) also implies that Cijkm = 
Ckmij- We subsume these relations as “pair permutation” symmetry. Physically, they 
correspond to the exchange of Cooper pairs. Moreover, as can be seen from reference 
tables for Clebsch-Gordan coefficients, we have Cijkm = to which we will 

refer as “pair sign-flip” symmetry. Physically, it is a consequence of the invariance of 
our system under reflection in the xy-plane. 

It thus suffices to look up (6.7) in a reference table for Clebsch-Gordan coefficients 
once for each member of a “pair permutation” and “pair sign-flip” equivalence class, 
ignoring those tuples (i, j, fc, m) which do not satisfy the selection rule i+j = k+m. The 
result is presented in Table 1. By counting the number of elements of each equivalence 
class, we find 

^fwave (^-2,..., ^2) 

2 

= ^{{'^ -t-+ 2 ^ 

V \m=—2 ) m=l,2 

- 21k ^0 + 23? ^^0 ^2^-2^ — 43? (^i'ip-i'il>2'ip-‘2^ \ , 


where r 


5^. Notice that this expression contains a second complete square: 


cGL 

^d-wave 



5c 

Idvr 



+2 


+ 2 ^ 2^-2 


( 6 . 8 ) 


To conclude Theorem 3.1, it remains to make the basis change to the real-valued 
spherical harmonics, i.e. to invert (3.2). On the level of the GL order parameters, this 
yields the SU{5) transformation 


■00 = V’o, 


^-2 


-V'l -I- fV’-i 

71 

^p2 - i'lp-2 

72 ’ 


•01 

02 = 


_ 01 -I- i0_l 

72 

02 -I- f0-2 

72 ■ 


(6.9) 
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6.3 Proof of Theorem 3.5 


The situation is as in three dimensions, only simpler. The dij GL coefficients are again 
diagonal by orthogonality and they come with a factor d defined in the same way as in 
Theorem 3.1 but with f 2 {p) replaced f 2 {p)/p since D = 2 (of course the definition of q 
has changed as well). For the Cijkm coefficients, instead of considering Clebsch-Gordan 
coefficients, it suffices to compute 

27r 

cos( 2(/9)^ sin(299)"^“^d(/? (6.10) 

for all 0 < A: < 4. Here, the GL coefficient c is defined in the same way as in Theorem 
3.1. We omit the details. 



6.4 Proof of Theorem 3.7 


We compute using the formulae (2.20) and (2.21) for the GL coefficients 

as in the previous section. We already computed most of the GL coefficients, namely 
all the ones that couple d-waves to d-waves. 

By orthonormality of the spherical harmonics, dij is still diagonal. For i, j 7 ^ s, d is 
as in (6.4). Notice however that d depends on q through / 2 . When i = j = s, we have 


to replace by Qs, which is conveniently described as multiplication by gs 


conclude that 


dij 


d( 2 «) [[ i = j = s, 

ddij otherwise. 


s± 

e 


. We 


with as defined in (3.23). 


We turn to the quartic GL coefficients Cijkm- Note that the “pair permutation” 
and “pair sign-flip” symmetries described in the proof of Theorem 3.1 still hold. In 
addition to the results listed in Table 1, we now have equivalence classes of Cijkm where 
some indices are equal to s. Since the corresponding dg carry zero momentum in the 
^-direction, the selection rule dictates that Cijkm can only be non-zero if the s replaces 
a 0 -index. 

We thus consider all equivalence classes of GL coefficients that can be obtained by 
replacing a 0 in Table 1 by s. We compute their values again via ( 6 . 6 ) (some follow 
immediately from the fact that Tg^ = l/\/4vr)- The results are presented in Table 2. 

Just as for dij, the are the result of a radial integration where for 

each index equal to s, is multiplied by a factor Qs- This yields the expressions 
(3.23) for Note that according to Table 2, CggsO = 0 and thus it is not 

necessary to define 

Armed with Table 2, it remains to count the number of GL coefficients in each 
equivalence class. After some algebra, we obtain 

^(s + d)-wave ^—2; • • • > V’2^ 

‘^d-wave(^—2) • • • ; V’ 2 ) T ^s-wavei^s) + “^coupling(^•S) ^—2; • • • > '02)- 
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i 

j 

k 

m 

^ijkm ' 287r 

s 

2 

0 

2 

-2^50^") 

s 

2 

s 

2 

7c(2+ 

s 

2 

1 

1 


s 

1 

0 

1 


s 

1 

s 

1 

7c(2+ 

s 

1 

-1 

2 

-^300^") 

s 

0 

0 

0 

2y5c(i*) 

s 

s 

0 

0 

7c(2+ 

s 

0 

s 

0 

7c(2+ 

s 

s 

s 

0 

0 

s 

s 

s 

s 

7c(4+ 

s 

0 

2 

-2 

-2^50^") 

s 

s 

2 

-2 

7c(2+ 

s 

0 

1 

-1 


s 

s 

1 

-1 

_7c(2+ 


Table 2: Equivalence classes of new Ginzburg-Landau coefficients in the mixed (s + (i)-wave 
case. are defined in (3.23). 


where 


^conp\ingi^s,i’-2, ■ ■ ■ ,^2) 

\/5cdd 


Ttt 


23ff 


I ^ ^ I'0m I 2 ^ ( I'0m I 

ym=0,itl m=±2 

— 2 ^ m3ff 

m=l,2 


050O0m0-r: 


+\/6 ^ 3 ^ 


cr=±l 


0<7 '0s'02o 


-23? 


0500 - 0-^020 


+ 


02^) 


^ r ^ 2 / ^ 

2|05p ^ |0imP + 3? 05 (^00 -2010-1 + 2020-2 


m=—2 


where <?^wave(0-2> • • • > 02) is given by (6.8) and 

^fwave(00 = ^ ((l0sp - ^-5) - 

with r5 = Statement (i) in Theorem 3.1, which gives the expression for 

+(i)_wave’ follows by transforming into the basis of real spherical harmonics 
via (6.9). 
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To prove (ii), we use the GL energy expressed in the basis of real spherical har¬ 
monics. Let e > 0 and take (V’-2, • • • £ A^d-wave, the set of minimizers of £d-wa.ve 

described by (3.12). Set ^ps = ecu with |a;| = 1 and note that 

+ d)-wave = inf Sd-Lve + e^pz] 

+ ^ 1 -— 1 +• 

\ TT 2tt J Att 

for some z G C, which is independent of e and w. Consider first the case that 
{ip- 2 , ■ ■ ■ ,'4>2) £ -Md-wave is such that z / 0. Then, we can choose u such that 
Repz] < 0 and we obtain (3.25) for sufficiently small e. Thus, suppose that z = 0, 
which is e.g. the case for (0, r/\/2, 0, ir/\/2,0) G -Md-wave- It is then clear that (3.25) 
holds iff ; oi' equivalently This proves (ii). 


For statement (hi), let V’s be a minimizer of S^-wave^ be. = Tg. Now let e > 0 
and let (ip- 2 , ■ ■ ■ P 2 ) have entries of the form iprn = eV’m with PP < 1- We have 


^{s + d)-wave (V’s) '0—2) 


GL 


= minTi 


+ 0{ep 


+ s^ 


,02) 


—d + 


42s), 


TT 




2^ C 
ml ' 


(2s) 

2p 




01 ^ (00)^ 


m=—2 


as e —)• 0. The real part is clearly minimal when we choose Aigp'P) = Aigps) + '7r/2 
for all m with ^ 0. This choice yields 


d)-wave(0^>0-2,---)02) 



+ 0{ep. 


When the term in parentheses is strictly negative, which is equivalent to < 

we see that £fs\d)-wave ^ ™™‘^sI0ave 1°^ Sufficiently small e. Vice-versa, when 
the term in parentheses is strictly positive, d)-wave ^ ™™‘^0wave lo^ small e > 0. 

To conclude statement (hi), it remains to consider the case when 

the 0(e^)-term vanishes. The leading correction is now given by the 0(e^)-term and 
by choosing ipm = 0 for m / 0, we find 


£GL 

^(5 + (i)-wave 


(05,0,0,00,0,0) 


min T 


GL 


+ £' 


;2v/5c(^0_„,2^r-w 


7tt 


I0ol 3f^[0s0o] + 0(e 


Letting Arg(V'o) = Arg(V’5 ) + tt shows that £fXd)-^^we (0s, 0, 0, 0o, 0,0) < minT^^^,^^^ 
in this case as well. This proves statement (iv). □ 


7 Proofs for part III 

The proof of Theorem 4.1 is based on three steps. 
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• In Lemma 7.1, we solve the eigenvalue problem for Kt + explicitly in each 
angular momentum sector 77;. The key result is the “eigenvalue condition” (7.3) 
which gives a formula for the eigenvalue (or energy) E in terms of the other 
parameters I, T and A. We will see that one can solve this for A and one obtains 
an integral formula which is monotone in E. Therefore, instead of showing that 
E is minimal for I = Iq, one can equivalently show that A is minimal for I = Iq. 

• In Lemma 7.2, we show how, by adapting the parameters in,T of the “weight 

function” p/Kt{p), one can conclude that Kt(p) positive, if one 

assumes that / is strictly positive on an interval. 

• By Theorem A.l, for any half-integer Bessel function of the hrst kind J^ig+ii 2 : 
there exists an open interval around its first maximum on which it is strictly 
larger than (the absolute value of) all other half-integer Bessel functions. 

The idea is then to use the eigenvalue condition (7.3) to rephrase the question whether 
some state in T-Li^ has lower energy than all states in 77; as the more tangible question 
whether the quantity 



is positive. By Theorem A.l there is an interval of p-values on which the integrand is 
positive and by Lemma 7.2 there are intervals of p- and T-values such that the entire 
integral is positive. 

7.1 Solving the eigenvalue problem 

For any radial V, we can block diagonalize Kt -|- F by using the orthogonal decom¬ 
position of L^(M^) into angular momentum sectors (3.5), namely L^(M^) = 0 ^q 77; 
with 77; dehned in (3.3). It is well-known [45] that the Fourier transform leaves each 
77; invariant. Consequently, if we have a € 77^(M^) satisfying the eigenvalue equation 

(Kt + V)a = Ea, (7.1) 

then we can decompose it as a = with a; G 77; mutually orthogonal. Taking the 

Fourier transform of (7.1) and using the fact that Vai G 77; since V is radial, we get 
from orthogonality 

KT(p)ai(p) + Vai(p) = Eai(p), 

for every I > 0 and a.e. p G Thus, we can study each component a; separately. 
When is the specihc radial potential (2.5), we can say even more. 

Lemma 7.1. Let Vx^r be as in (2.5) and let I be a non-negative integer. We write 
for the Bessel function of the first kind of order I + 1/2. Let E < 2T if p > 0 

E < tanh(|!l|/(2T)) < 0. Then 

ker (Kt + Vx,R - .E) n 77; / 0 (7.2) 
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is equivalent to the “eigenvalue condition 


1 = A 


pR 


-j‘^MpR)dp. 


/o KT{p)-E"i+k 

Moreover, if (7.2) holds, then ker {Kt + Vx^r — E) = span{/9;} (g) Si with 

roo Ji^iirp)Ji^i(Rp) 


(7.3) 


roo ^ 

Pl{r) = r~^l‘^ / p- 
JO 


Kt{j>) - E 


-dp. 


(7.4) 


Since \f^l_^_l{p)\ < Cp the numerator in (7.3) and (7.4) poses no threat for 
convergence of the integral. 


Proof. By the definition of Ri, we have 

i 

«z(x) = ^ cxi^m{r)Y{^{'d,(f), yL={r,d,p). 

m=—l 

We suppose a; satisfies {Kt + Vx^r) og = Eai. Recall that the Fourier transform not 
only leaves each Ri invariant, it also reduces to the Fourier-Bessel transform Ei on it 
[45]. That is, a function of the form /(x) = g{r)Yf^{'d, (p) has Fourier transform given 

f{p) {Y'ig){p)Yi^{P,ip), p = {p,P,ip), (7.5) 

where the Fourier-Bessel transform reads 

roo 

J^ig{p)= s^/‘^p~^/‘^Ji^i{sp)g{s)ds. (7.6) 

Jo 2 

We apply the Fourier transform to the eigenvalue equation. By (7.5) and orthogonality 
of the spherical harmonics, 


{Kt{p) — E)Eiai^rn{p) + Ei{V\^Rai^m){p) = 0 


for all m and a.e. p G M^. 
therefore 


The assumption on E is such that Kt{p) 


ElCXl^m{p) — 


JEi{V\fiai^m){p) 

Kt{p) - E 


E > d and 
(7.7) 


So far we only used that the potential is radial. Since Va,r = —^^{\ ' I ~ 7?), 


-J^l{yX,RC^l,m){p) = -ai,m{R){ElVx,R){p) 

= Aaz,™(R)R3/2p-^/2j-^+i(Rp). 

Plugging this back into (7.7), we find the following explicit expression for the solution 
to the eigenvalue problem; 


(p) 


Aaz,^(R)R3/2p-i/2 


Ji^{Rp) 

Kt{p) - E 


(7.8) 
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Now we apply ^ which, by unitarity of the Fourier transform, has the operator 
kernel when evaluated at r > 0. For all r > 0, we have 


roo 

ai,mir) = / p 

Jo 


~ Ji+Urp)Ji+i{Rp) 


Kt{p) - E 


dp 


Note that we may assume that for some m, ai^rn{R) / 0, since otherwise a; = 0. 
Evaluating the above expression for that particular m at r = R gives (7.3). We write 
Oii^m{R) = and absorb ci^m into the angular part Si to get (7.4). Clearly 

the argument works in reverse, proving the claimed equivalence. □ 


7.2 Choosing /i and T 

From now on, let /U > 0. The following lemma concerns the quantity 



V 

Kt{p) 


f{p)dp. 


Suppose we know that f > £ on some interval I, while / may be negative outside of I. 
Our goal in this section is to choose the right values of p and T such that the above 
integral is then also positive. 

The basic idea is to view p/Kt{p) as a weight function which is centered at the 
point p = where it takes a value proportional to T~^. By making T small enough, 
we can ensure that the neighborhood of the point p = dominates in the above 
integral. By choosing G I and T sufficiently small, the integral will pick up mostly 
points where / is positive and will therefore yield a positive value itself. This is spelled 
out in the following lemma. 

We will eventually apply this lemma with / = ~ ■^i+i /2 positivity of 

the above integral will translate via (7.3) to the statement that the angular momentum 
sector Ffzg has lower energy than Ri. 


Lemma 7.2. Let f : M+ —)• M 6e a continuous function satisfying \f{p)\ < Cf{l + 
p2^-i/2 some Cf > 0. Suppose there exists £ > 0 and an interval {a,b) such that 
f > £ on {a,b). Then, 

(i) for every (5 > 0 small enough, there exists T* > 0 and an interval I such that for 
every /U G / and T G (0,r*), 

u2 2 

(a) letting S := , one can choose 


I := + s,b^-6), 



2Cf (yi + 262 + \ 

I ■ 


(7.10) 


44 



Proof. Let /U G (a^ + 5, 6^ — 5). Since > 0 and | tanh | < 1, we can estimate 

ri^^dp>-c,/ j p.ii) 

Jo Kt{p) ■' y[0,a)u(b,oo) (1 Ja Kt{p) 

In the first integral, we estimate pointwise 

\P^ - P\~^ < S~^{X{p<2b} + ‘Jp~^X{p>2b}) 

with XA denoting the characteristic function of a set A. This gives 

P 


-C 


f 


[0,a)U(b,oo) (1 +p2)^/2|p2 _ ^1 


dp > -Cf6-^ (\/l + 262 + (7.12) 


In the second integral, we change variables and use p G {of+6, b‘^ — 5) with tanh(u)/u > 
0 to get 


P 


Kripf^ 2 


1 /” 2 t tanh(u) 


2T 


du > 
u I A 


2 T tanh(u) 


u 


du > log ( — 


where in the last step we also used that tanhx >1/2 for x > 1. Combining everything, 
we get 


Pfip) 


Jo Kt{p) 

The claim follows from some algebra. 


dp > -Cf5 ^ ( \/l + 262 + _ I + £log ( — 


(5 


2T 


(7.13) 

□ 


7.3 Proof of Theorem 4.1 


Proof of (%). By rescaling the parameters p, X and T, we may assume that R = 1. We 
fix a non-negative integer Iq and invoke Theorem A.l to get e > 0 and an interval (a, 6) 
on which jf^j^Yi 2 ~ ■^i+i /2 ^ ^ ^ / ^ 0 - Then we apply Lemma 7.2 to 


/ •— ^«o+l /2 ■^/+l/ 2 > 

which satisfies 

|/(p)| 

and so Cj = 2^/2 Lemma 7.2. To prove (7.14), ones uses statement (ii) in Lemma 
A.5 to get Jy{p) < MJ;{p) < 1 for all v. Together with \Jy\ < 1 from (9.1.60) in [1], 
this implies \f{p)\ < min{l,p“^} and hence (7.14). 

Note that T* and I defined in Lemma 7.2 (ii) work for all I ^ lo, because they 
depend on / only through (a, 6), which is uniform in / by Theorem A.l, and through 
Cf = 2^/2 _ Hence, Lemma 7.2 provides T* > 0 and an interval I such that for all p £ I, 
all T < and all I ^ Iq we have 


r P 
Jo Kt{p) 


{•^l+i/2{p) - dlj+i/2{p)) dp > 0 


(7.15) 
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For every non-negative integer I, we define the function 

‘ (7.16) 

which is chosen such that A satisfies the eigenvalue condition (7.3) with E = 0. We 
write 

Ei{T,iJ,X) := inf spec (i^T + Va,i) 

With these definitions, Lemma 7.1 says 

Ei{T,^i,Xi{T,fi)) = 0 (7.17) 

At the heart of our proof is the following monotonicity argument. For all n £ I, all 
T <T^ and all I Iq, we have 

0 = Ei{T,^i,Xi{T,fi)) < Ei{T,fi,Xi^{T,iJ,)), (7.18) 

where the inequality holds by the variational principle applied to the operator {Kx + 
and the observation that (7.15) is equivalent to Xig{T,fi) < Xi{T,fi). (The 
inequality is strict because {a, VApa) = —Xa{Kf' is either strictly monotone decreasing 
in A or identically zero and in the latter case the energy has to be at least 2T.) 

This would already prove (4.2) and (4.3) under the condition that one fixes T <T^ 
and determines A through (7.15). We find it physically more appealing to fix A small 
enough and determine T instead. To this end, we observe that T i—)■ A^p (T, //) is 
monotone increasing, because T i—)■ Kt{p) is monotone increasing for every p > 0. 
Therefore, for every p £ I, we have the monotone increasing inverse function 

{Oj p)) —> (0)7"*) 

A I—)■ r(A, p) 

satisfying Xig(T{X, p), p) = X. To remove the ^u-dependence from the maximal value 
for A, we set 

A* ;= minA(r*,//) (7.19) 

l-t^I 

and note that A* > 0 since the integral in (7.16) is continuous in p by dominated 
convergence. For A < A*, (7.17) and (7.18) become 

Eio{TiX,p),p,X) = 0, Ei{T{X,p),p,X) > 0, VI ^ Iq. 

This proves that for all p £ I and all A < A*, there exists Tq < T* (namely Tq := 
T{X,p)) such that (4.2) holds (modulo restoring the R parameter). Moreover, (4.3) is 
a direct consequence of the explicit characterization of ker(iF7’c + in Lemma 7.1. 
Finally, (4.4) follows via the variational principle from the observation that T i—)■ Kt{p) 
is strictly increasing for all p > 0 and so T i—)■ Eiq{T, p, A) is strictly increasing as well, 
as long as it stays below 2T. □ 
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Proof of (ii). Consider the function 


Claim; There exists T** > 0 such that for allO < T < T** there exists > 0 such that 
hrihr) = 0. Moreover, y/JTr —> ^i /2 as T —)■ 0, where Zxj^ = min{z > 0 : = 

The claim follows essentially from the intermediate value theorem. Before we give 
the details, we explain how one may conclude statement (ii) from the claim. Let 
0 < T < T**. By definition (7.16), 5t{tt) = 0 implies \q{T,^t) = ^ 2 {T,ij,t)- By 
Lemma 7.1 and using the notation (7.17), 

Eo{T, IJ,T, Xo{T, ht)) = -E'2(T,/ir, Ao(r, ^r)) = 0. (7.20) 

This implies C in (4.7) according to Lemma 7.1. Equation (4.8) follows by the same 
monotonicity argument as in the proof of statement (i) above. 

In order to prove (4.6) with the choices /x = /xy and A = Ao(T, /xt) and the remaining 
D in (4.7), we shall show that there exists T* G (0,T**] such that for all 0 < T < T*, 

£'/(r,/XT, Ao(T,/X t)) > 0, V^>4, Hseven. (7-21) 

By Theorem A.l (ii) (with lo = 1) and Lemma A.6, there exists an open interval 
containing zi /2 such that 

■^ 5 % “ 77 ^+i/2 > interval. 

I even 

As in part (i), Lemma 7.2 provides T** > 0 and an interval I' containing such that 
for all /X G all T < T** and all even / > 4 we have 

[ 1 ^) ~ ^ 

Since the second part of the claim gives /xt —)• Zy^ as T —)■ 0, we may assume, after 
decreasing T** to T* if necessary, that ixt G I' for all 0 < T < T*. Therefore (7.22) 
implies that Ao(T, /xt) = A 2 (T, /xt) < for all T < T,, and all even I > 4. By 

the same variational argument as in (7.18), this implies (7.21). 

We now prove the claim. The reader may find it helpful to consider Figure 1. 
Since /x i—)■ Kt{p) is continuous for every p, p ^ St is also continuous by dominated 
convergence. Let xi {I = 0,2) denote the first maximum of 77;_|_i/2- It is well-known 
that xq < X2 [37] and that > J^i2{xo) and (which is also 

a very special case of our Theorem A.l (i)). By continuity these inequalities hold also 
in neighborhoods of xq and X 2 . Therefore Lemma 7.2 provides open intervals Ii C M+ 
{I = 0,2), containing xi, and a T** > 0 such that for all T < T**, we have (5 t > 0 on 
Iq and 5 t < 0 on / 2 - By the intermediate value theorem, for any T < T** there is a 
/XT £ [sup Iq) iaf 72 ] with St{pt) = 0. This proves the first part of the claim. 
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We are left with showing that ^/JIT —^ 2:^/2 as T —)■ 0. Since £ [sup/o; inf/ 2 ] 
is bounded, it has a limit point as T —)■ 0. We argue by contradiction and assume 
that there is a limit point z different from Zi/ 2 - By Lemma A. 6 , Zi /2 is also the 
position of the first critical point of By the interlacing properties of the zeros 

of Bessel functions and their derivatives, see e.g. [37], Zi /2 £ {xo,X 2 ) and there is no 
other point z G (xo,X 2 ) at which = ^‘ii 2 ^^)- Therefore J‘li 2 ~ either 

strictly positive or strictly negative at z and, by continuity, also in an open interval 
containing z. Lemma 7.2 provides an open interval 1 containing z^ and a T > 0 such 
that is either strictly positive or strictly negative for all T < T and /j, G I. Since 

z is a limit point of there is a sequence Tm 0 with —)■ z^. In particular, 

UTm ^ L and Tm < T for all sufficiently large m. Thus, is either strictly 

positive or strictly negative for all sufficiently large m. This, however, contradicts the 
construction of jiT, according to which (5r(MT) = 0 for all T < T**. Thus, we have 
shown that ^/JIT —^ 21 / 2 - Tl 

A Properties of Bessel functions 

While one might expect the following fact about Bessel functions to be known, it 
appears to be new: 

At its first maximum, a half-integer Bessel function is strictly larger than (the 
absolute value of) all other half-integer Bessel functions. 

The precise statement is in Theorem A.l (i) below. It extends to families of Bessel 
functions {Ju+k}k&’i+ with n G [0,1], in particular to the family of integer Bessel 
functions. We acknowledge a helpful discussion on mathoverflow.net [36] that led to 
Lemma A.5. 

Let Iq be a non-negative integer. We recall that the Bessel function 77 ;g+i /2 (of 
the first kind, of order Iq -|- 1 / 2 ) vanishes at the origin and then increases to its first 
maximum, whose location we denote as usual by \ - The following theorem says 

that at 1 ’ ■^io+i /2 strictly larger than any other with I a non-negative 

integer different from Iq. 

Theorem A.l. Let Z+ denote the set of non-negative integers and let Iq G Z+. Recall 
that j[^j^]^j 2 1 denotes the position of the first maximum of 12 - 

(i) There exist e > 0 and an open interval I containing 1 that 

^Zo+ 1 / 2 - sup J' 4 i/ 2 >e on I. (A.l) 

zez+\{Zo} 

(ii) Iflo > I, then Jia-i/ 2 Ul+i/ 2 ,i) = ^Zo+ 3 / 2 (/((,+i/ 2 ,i) t/iere exist e' > 0 and 
an open interval I' containing f[^^Yi 2 1 ^^oh that 

min{J; 2 _^/ 2 >^Zo+ 3 / 2 } - sup jf+^i 2 > ^ on l'. (A.2) 

l—lo odd 

Remark A.2. Statement (i) is the key result and implies Theorem 4.1 (i). Statement 
(ii) is used to prove Theorem 4.1 (ii). 
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Figure 1: A plot of the squared Bessel functions Jlj 2 ^ ’ •^ 2 ^ 2 - Observe that in 

an open interval around its maximum, each function is the largest one among all the shown 
ones (in particular it is the largest among all the according to Lemma A.3). 


The proof of (i) in Theorem A.l is split into three Lemmata, each treating one of 
the following three regimes of 1: 

jC> : = {? G ; I > Zq} ) 

^ '^+ ■ ^ <^0, ji+1/2,1 > Jzo+1/2,1} ’ 

£< : = G Z+ : I < lo, ji+i/ 2,1 < jzo+ 1 / 2 , 1 } • 

Here, as usual, ji+i/ 2,1 denotes the first positive zero of t7/+i/2- The most cumbersome 

regime is £<g. The proof there is based on a combination of some hands-on elementary 
estimates and bounds on the zeros of Bessel functions and their derivatives, which we 
could not find in the usual reference books [1],[50]. The first regime £> is the easiest 

Lemma A.3. There exist ei > 0 and an open interval Ii containing 1 such that 

^Zo+ 1/2 - sup ^41/2 ^ on h (A.3) 

1>Iq 

Proof. According to [31], the function 

n max|X(2/)| 

y 


is strictly decreasing. Therefore 


•“ 2 (^L+i/2(-?Zo+i/2,i) 


max^. 


Zo+3/2 
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is strictly positive. By continuity, there exists an open interval Ii containing i 

such that for all x G Ii, 


For X G Ii, we have 

^Zo+i/ 2 (®) - sup J’4i/2(2;) > -£i + ^io+i /2 (4+1/2,! ) “ sup max /7-4i/2 (?/) 

l>lo l>lo y 

— “^1 + ■^41/2(4+1/2,1) “ “■^^■^43/2(1/) 

= ei- □ 


Lemma A.4. There exist £2 > 0 and an open interval I 2 containing \ such that 


4+1/2 - sup Ji +^/2 > £2 on h. 


Proof. Since the supremum of finitely many continuous functions is itself continuous, 
it suffices to prove 4+i/2(4+i/2,i) > ^z4/2(4+i/2,i) ^ ^ 4' define 

the sequence {az}ze/;< by 

^Z+1/2(JZo+1/2,i) = ®z4+1/2(/Zo+1/2,i)- (^-4) 

With this definition, the recurrence relation for Bessel functions from (9.1.27) in [1] 
appears in the form of a second-order difference equation 


I + 1/2 

oz-i = 2- ai — az+i 

Xq 


(A.5) 


with initial conditions ai^ = 1 and ai^^i = (Zq + l/ 2 )/j 4 i /2 i' well-known that 
the latter quantity is strictly less than one, see eq. (3) on p. 486 of [50]. Moreover, 
az > 0 for all I G C<, because jz+i/ 2 ,i > 4+1/2 i Bessel functions are positive 

before they first become zero. An easy induction lets us conclude from (A.5) that 
ai < az+i < 1 for all I G C<. In particular, ai < ai^-i = (Zq + l/2)/4+i/2 i ^ 
Recalling the definition (A.4) of a;, this proves the claim. □ 


We finally come to the regime £<g. As a tool, we will use the “modulus” function 
defined by 

M, := 44 + 4 , 

where yu is the Bessel function of the second kind. The first two statements of the 
following Lemma are known facts about the modulus function. Statement (iii) is the 
key result to derive (iv). 


Lemma A.5. (i) The map u i-+ My{x) is strictly increasing for all x > 0. 


(ii) For all x > v, 


M/(x) < 


2 1 
"Zl \/x^ — 
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(Hi) If lo > 11, there exists li < Iq such that we have both, 

(Hi.a) ^Jq+i/20’zq+1/2,i) > ^h+l/2^ilo+l/2,l) 

(iH.h) jzi+l/2,1 > 4+1/2,1 

(iv) There exist 63 > 0 and an open interval I 3 containing 4+1/2 1 that 

Jl+i /2 - sup > ^3 on I 3 . 

The intuition why such li as in (iii) should exist is based on a heuristic argument of 
which we learned through [36], involving asymptotic formulae for the relevant expres¬ 
sion. To turn this into a rigorous proof, we need to replace the asymptotics by bounds 
that hold for all Iq (or at least for alHo ^ H)- [18] contains results which are sufficient 
for our purposes when combined with a number of elementary estimates. 


Proof. Statement (i) is a direct consequence of Nicholson’ formula, see p. 444 in [50], 
and the fact that Kq > 0. Statement (ii) is formula (1) on p. 447 of [50]. 

We come to statement (iii). For convenience, we write m = I + 1/2, so in partic¬ 
ular mo = ^0 + 1/2. We also abbreviate xq = 4-i-i/2i- basic idea (inspired by 
asymptotics) is to choose 

1/3 

mi = mo — cmQ 


with c small enough to have (iii.a) hold but large enough to have (iii.b) hold. By (i), 
(iii.a) is implied by 


2 1 

TT yjxl-ml 




(A.6) 


By [18], we have the lower bound 


xo > mo exp 





i.oemg^/^) 


(A.7) 


for all mo > 11.5. Here, a'l is the absolute value of the first zero of the derivative of 
the Airy function, with a numerical value of about 1.018793. From mo > 11.5, we 
can conclude that the argument of the exponential in (A.7) is greater than O.Gmg ' . 
Thus, by the elementary estimate > 1 + y, (A.7) implies the more manageable lower 
bound 

Xo > mo -|- 0.6mg'^^ 

1 /3 

Setting mi = mo — cmg' with c to be determined and using the above bound on xo, 
as well as mo > 11.5, we see that (A.6) is implied by 


2 1 
vr Vl.26-h 2c - 0.19c2 



max| Jmo(x) 

X 


2 


(A.8) 


According to [31], u 1—)■ ma,^x\IIv{x)\ is an increasing function and so we can 

estimate the right-hand side in (A.8) from below by Jyixf' for any 1/2 < 

V < mo. Unfortunately, the numerical value one obtains for the “worst case” v = 1)2 
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is not good enough to also get (iii.b). Instead, we assume that c < 1 and use mo > 11-5 
to get mo — crug > 8.5 and so 

(ml^^max\Jmo{x)\^ > (^(8.5)^/^ max | J8.5(aj)|) >0.42 

where the last inequality can be read off from a plot, for example. Therefore, (A.8) 
holds if we can find c < 1 that satisfies 


2 1 
vrVl.26 + 2c-0.19c2 


< 0.42 


(A.9) 


and it is easily seen that this holds for c G [0.5,1]. 

Now, we want to ensure that c is also small enough to have (iii.b) hold, i.e. jmi > xo- 
To this end, we invoke two more facts: 

• the upper bound 

Xo < mo + O.SOm^^. (A. 10) 

This is a consequence of the bound 


Xo < mo exp 


(2 


from [18], where again ~ 1.018793, by noting that mo > 11.5 implies that the 
argument of the exponential, call it y, satisfies y < 1.59. On [0,1.59], we can 
estimate exp(y) < 1 + 1.09y, as one can verify e.g. by plotting and this yields 
(A.IO). 

• the lower bound 

jmi > mi + (A.ll) 

which we obtained from the optimal lower bound proved in [38] by rounding 
down. This is better than the bound one can derive from a corresponding result 
of [18] as we did above. 

From (A.IO) and (A.ll), we see that jmi > xo will follow from 

(mo — cmy^) + 1.85(mo — cmg'^^)^/^ > mo + 0.89mg'^^, (A.12) 

Since c < 1 and mo > 11.5, we have 1 — cmg > 0.8 and so (A.12) is implied by 

c < (0.8)^ * 1.85 - 0.89 = 0.827. 


So any choice of c G [0.5,0.8] will ensure that (iii.a) and (iii.b) hold. 

We prove statement (iv). By continuity, it suffices to prove 77;g_|_i/2(*o) > J^{xo) 
for all I G £< (which we recall means I < lo with < xo)- Assume first that 

lo > 11. Choosing li as in statement (iii), (iii.a) states 

■^io+i/2(^o) ^ ^h+i/2i^o) (A.13) 
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and (iii.b) implies that li G C<. By the monotonicity oi v ^ j^, it holds that I G 
implies I < li. Thus, the definition of and statement (i) imply 

-^Z+1/2 — -^Z^+1/2- (A-14) 

Together with (A. 13), this implies (iv) for Iq > 11. Since for Iq = 0,1 there are no 
I <C lo, we may assume Iq > 2. For 2 < Iq < 10, one can then check by hand that 
(A. 13) holds with the choice li = Iq — 2. Since Iq — 1 G C<, we get that I G £< implies 
I < li and so (A.14) applies for all such 1. □ 

Lemma A.6. For any positive integer I, 

min{z > 0 ; J^_ii 2 {z) = J/+ 3 / 2 (^)} = fi+i/ 2 ,i (A.15) 

and Ji_ii 2 ,Ji+ii 2 ,Ji+ 3 l 2 are positive on (0, 

Proof. We recall the recurrence relation from (9.1.27) in [1], which says that for all 
i',z > 0, 

Jiy-i{z) - Ju+i{z) = 2jf{z). 

Applying this with n = / + 1/2, z = obtain Jz-i/ 2 (j)'+i/ 2 ,i) = Ji+ 3 / 2 {j'i+ii 2 p) 

and hence < in (A.15). Notice that by the interlacing properties of zeros and extrema 
of Bessel functions, see e.g. [37], i positive zeros of 

77i-i/2) 77 ;+i/ 2, 77;+3/2- Since Bessel functions are positive before they reach their first 
positive zero, we conclude that 77;_i/2,77;+i/2, 77)+3/2 are positive on (0, il- 
particular, 77;_i/2, 77/+3/2 are positive at the left side of (A.15), call it z;, and so we 
can take square roots to get Ji-i/ 2 {zi) = Ji+ 3 i 2 {zi). By the recurrence relation from 
above, 77’/^_3^2('^0 ~ ^ implying zi > j[_|_i /2 i’ claimed. □ 

It remains to give the 

Proof of Theorem A.l. Statement (i) is a direct consequence of Lemmata A.3 to A.5. 
For statement (ii) we first observe that for any positive integer I, 

■^Z-l/2 > ■^Z+3/2> (0) Jz+l/2,l)- (A.16) 

In fact, by standard asymptotics, this inequality holds near zero and, according to 
Lemma A.6, J/+ 1/2 1 point of intersection of Ji_yi 2 ^^*4 ■^z+ 3 / 2 - Therefore 

the inequality holds on all of (0, j[_|_^y 2 i); claimed. 

We now use the fact that j'i_^_ij 2 1 is increasing in I [37]. Choose I' to be an open 
interval containing j[(j_|_i /2 1 whose closure is contained in (0, i)- Then by (A.16) 

(with I = Iq + 2) and continuity there is an e' > 0 such that 

■i^Zo+3/2 - '^^0+712 + ■ 

Applying (A. 16) successively with / = /q + 4, /q + 6,..., we conclude that 

^Zo+3/2> sup J'4i/2 + e' on/', 

;>;o+3 
l—lo odd 
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which is one part of the claim. Finally, we want to prove the same inequality with 
■^io-i/2 (with possibly smaller e' and I'). Clearly, (A.16) implies that 

this is true on I' n (0, il- Now use continuity to find <5 > 0 such that > 

■^io+3/2 “ e /2 on [j)o+i/2,i’il)o+i/2,i 

/t2 \ /T’2 I // 7"// 

^Iq-\I2 ^ ^Zo+7/2 + ^ 

with e" = eY2 and 1” = /' n (0, i + before, (A.16) now implies the 

inequality in part (ii). This completes the proof. □ 
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